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ABSTRACT 

The  scattering  of  a  scalar  plane  wave  by  a  totally  reflecting  sphere 
(hard-core  potential)  at  high  frequencies  is  treated  by  a  modified  Watson 
transformation.   The  behavior  of  the  solution  both  in  the  near  and  far  regions 
of  space  is  discussed,  as  well  as  the  accuracy  and  domain  of  applicability 
of  the  WKB  approximation  and  classical  diffraction  theory.   It  is  shown 
that  different  transformations  are  required  in  the  forward  and  backward 
half -spaces,  and  corresponding  integral  representations  for  the  primary 
wave  are  derived.   The  transformations  are  rigorously  proved  and  the  con- 
vergence of  the  residue  series  is  discussed.   In  the  shadow  region,  the 
physical  interpretation  of  the  complex  angular  momentum  poles  in  terms  of 
surface  waves  is  in  agreement  with  Keller's  geometrical  theory  of  diffraction. 
In  the  lit  region,  sufficiently  far  from  the  shadow  boundary,  the  WKB  expansion 
for  the  wave  function  is  confirmed  up  to  the  second  order.   On  the  surface 
of  the  sphere,  Kirchhoff's  approximation  is  accurate,  except  in  the  penumbra 
region,  where  the  behavior  is  described  by  Fock's  function.   The  diffraction 
effects  in  the  neighborhood  of  the  shadow  boundary  are  investigated  and  the 
corrections  to  classical  diffraction  theory  are  obtained.   The  shift  of  the 
shadow  boundary  is  evaluated.   The  expression  for  the  wave  function  in  the 
Fresnel-Lommel  region  is  derived  and  applied  to  the  discussion  of  the  Poisson 
spot  and  the  behavior  near  the  axis.   The  total  scattering  amplitude  is 
evaluated  for  all  angles,  including  the  neighborhood  of  the  forward  and  back- 
ward directions.   The  corrections  to  the  forward  diffraction  peak  and  the 


transition  to  the  region  of  geometrical  reflection  are  discussed.   The 
modified  Watson  transformation  is  also  applied  directly  to  the  scattering 
amplitude.   The  connection  between  representations  valid  in  different  regions 
is  established. 


I.   Introduction 

The  problem  to  be  considered  is  the  scattering  of  a  scalar  plane 

wave  by  a  totally  reflecting  sphere  at  high  frequencies.   This  means  not 

1/3 
only  that  ka  »  1,  but  also  that  (ka)  '  »  1,  where  k  is  the  wave  number 

and  a  is  the  radius  of  the  sphere.   It  will  be  attempted  to  make  the  treat- 
ment both  rigorous  and  comprehensive,  including  a  discussion  of  the  near 
and  far  regions  and  a  comparison  of  the  results  with  the  WKB  approximation 
and  classical  diffraction  theory. 

The  wave  function  can  be  interpreted  either  as  the  velocity  potential 
of  sound  waves,  corresponding  to  an  acoustically  soft  sphere,  or  as  the 
Schrbdinger  wave  function  in  nonrelativistic  quantum  mechanics,  in  which 
case  it  corresponds  to  a  hard-core  potential.   In  both  cases,  the  boundary 
condition  is  the  vanishing  of  the  wave  function  on  the  surface  of  the  sphere. 
There  is  no  difficulty  in  extending  the  treatment  to  a  vector  wave  field, 
so  as  to  represent  electromagnetic  scattering  from  a  perfectly  conducting 
sphere. 

This  is  perhaps  the  simplest  problem  involving  a  finite  size  scatterer, 
and  its  exact  solution  in  the  form  of  a  partial -wave  series  has  been  known 
for  a  long  time.   It  is  also  well  known  that  this  form  of  the  solution 
becomes  useless  at  high  frequencies,  because  of  the  large  number  of  terms 
one  would  have  to  keep  in  order  to  get  a  good  approximation. 

The  way  out  of  this  difficulty  was  proposed  by  Watson(_l)  ,  who  trans- 
formed the  partial -wave  series  into  a  "residue  series",  which  is  rapidly 


convergent  at  high  frequencies.   Several  applications  of  this  transformation 
to  the  theory  of  radio  wave  propagation  around  the  Earth  were  made  by 
Van  der  Pol  and  Bremmer  (2) . 

Watson's  procedure  applies  only  to  the  shadow  region  behind  the  sphere, 
and  not  to  the  lit  region.   An  extension  of  the  treatment  to  the  lit  region 
was  made  by  Fock  (_3)  and  later  reformulated  by  Franz  (|+,^). 

As  will  be  shown  in  Section  VI,   this  treatment   is  also  incomplete, 
for  it  does  not  apply  to  the  backward  half -space.   A  certain  amount  of 
confusion  seems  to  exist  in  the  literature  concerning  the  application  of 
Watson's  transformation  in  the  lit  region. 

Recently,  interest  in  Watson's  transformation  has  been  renewed,  in 
connection  with  Regge ' s  work  on  complex  angular  momentum  in  potential 
scattering  and  its  applications  to  elementary-particle  physics  (6)  . 

In  view  of  this,  as  well  as  of  its  intrinsic  interest,  a  reexamination 
of  the  totally  reflecting  sphere  problem  seems  warranted.   Since  this  is  the 
simplest  problem  of  its  kind,  it  should  also  serve  as  a  model  for  the  exten- 
sion to  more  complicated  situations,  such  as  different  refractive  indices 
and  different  shapes.   The  case  of  a  transparent  sphere,  which  has  numerous 
applications  to  optics  and  nuclear  physics,  will  be  treated  in  a  subsequent 
paper. 

In  the  present  work,  Wa.tson's  transformation  will  be  reformulated  in 
such  a  way  that  it  becomes  applicable  both  in  the  forward  and  in  the  backward 
half -space,  and  in  particular  in  the  neighborhood  of  the  forward  and  backward 
directions.   For  this  purpose,  new  integral  representations  of  the  primary 
wave  will  be  derived  in  Section  II. 


A  rigorous  proof  of  Watson's  transformation  and  the  convergence  of 
the  residue  series  will  be  given  in  Sections  III  and  IV.   The  physical 
interpretation  of  the  results  in  the  shadow  region  will  be  discussed  in 
Section  V.   In  Section  VI,  the  nearby  lit  region,  excluding  the  neighbor- 
hood of  the  shadow  boundary,  will  be  considered.   It  will  be  shown  that  the 
solution  in  this  region  agrees  with  the  WKB  approximation  up  to  the  second 
order. 

The  diffraction  effects  arising  in  the  transition  between  lit  and 
shadow  regions  will  be  examined  in  Section  VII,  where  the  results  derived 
from  the  exact  solution  will  be  compared  with  classical  diffraction  theory. 
The  shift  of  the  shadow  boundary  will  also  be  investigated. 

At  greater  distances  from  the  sphere,  but  still  in  the  near  region, 
the  wave  function  resembles  the  classical  solution  to  the  problem  of 
diffraction  by  a  circular  disc,  as  will  be  shown  in  Section  VIII.   The  well- 
known  Poisson  spot  effect,  as  well  as  the  behavior  near  the  axis,  will  also 
be  studied.   Finally,  in  Section  IX,  expressions  for  the  scattering  amplitude 
in  all  directions  will  be  derived. 

The  main  results  and  conclusions  derived  from  the  present  treatment 
will  be  summarized  in  Section  X.  Some  of  them  are  not  new,  but  have  been 
included  here  for  completeness.   The  basic  mathematical  tools  employed  will 
be  presented  in  Appendices  A  to  F. 


II.   Watson's  Transformation 
A.   The  Total  Wave  Function 

Let  the  incident  plane  wave  be  given  by 

00 

t^(r,0)  =  exp(ikr  cos  0)  =  V  (2i  +  1)  i"^  j^(kr)  P^(  cos  9),  (2.1) 

where  J   is  the  spherical  Bessel  fianction  of  order  i  and  P. (cos  9)  is  the 
i       Legendre  polynomial. 

The  solution  of  our  problem  is  given  by  the  well-known  partial -wave 
expansion 


*(r,0)  =  \    E  (2i+l)i^ 
^=0 


hf^(p)  +  S/p)h^^^p) 


P^(cos  0),(2.2) 


where  the  S-function  S  ( p)  is  determined  by  the  boundary  condition  t(a,0)  =  0: 

S/p)  =  -\!>p{^)l\i^p{^),  (2.3) 

and  we  have  introduced  the  notations 

p  =  kr,    p  =  ka  (2.U) 

It  is  well  known  that  at  high  frequencies  one  can  associate  with  the 
i       partial  wave  an  "impact  parameter" 

p^  =  (i  +  ^)/k  (2.5) 

and  all  partial  waves  with  P„  ~  a  are  strongly  affected  by  the  scatterer. 
Thus,  the  number  of  terms  one  has  to  keep  in  the  partial -wave  expansion  is 
of  the  order  of p »  1,  so  that  (2.2)  becomes  useless  at  high  frequencies. 


Watson's  transformation  is  based  upon  the  following  formula: 

tfU^i)-'^^fMZi'l"^      ^'  (2.6) 


£  =  0 


where  C  is  the  contour  shown  in  Fig.  1.   This  formula  can  easily  be  checked 
by  taking  the  residues  of  the  integrand  at  the  physical  (half -integral)  values 
of  X. 

Clearly,  there  is  a  large  degree  of  arbitrariness  in  (2.6).   For 
instance,  the  factor  (cos  n\)        might  have  been  replaced  by  any  analytic 
function  having  poles  with  the  same  residues  as  this  factor  at  the  physical 
points.   The  only  restrictions  to  which  the  "interpolating"  function  f(x) 
is  subject  are  that  it  must  reproduce  t{£   +2)  at  the  physical  points  and 
that  it  must  be  regular  in  a  neighborhood  of  the  real  axis,  so  that  the 
integral  can  be  computed  by  residues  in  the  indicated  manner.   There  is 
usually  a  wide  class  of  functions  satisfying  these  conditions. 

The  choice  of  f(X)  is  dictated  in  practice  by  the  requirement  of 
appropriate  behavior  at  infinity  in  the  X  plane,  since  the  next  step  in 
Watson's  transformation  will  consist  in  the  deformation  of  contour  C  away 
from  the  real  axis.   In  Field  Theory,  this  leads  to  a  unique  continuation, 
with  the  help  of  Carlson's  theorem  (6). 

This  result  does  not  apply  in  the  present  case.   However,  only  two 
alternative  choices  of  f(x)  will  be  required.   The  first  one  leads  to 

nr,9)    =y  f(x,p,p)  P^_:  (cos  9)    e'^^^/^  ^^  ,  (2.7) 

where 
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with 


f(X,P,p)   =  |(.f-) 


1/2 


-ijt/i+ 


m^'^^"'^' 


;(^,3,p) 


H^-O) 


g(x,p,p)  =  h[^)(p)h[2)(p)   -h^2)^p^^i)(p^ 


(2.8) 


(2.9) 


In  this  expression,  R     i_     is  Legendre '  s  function  of  the  first  kind  and  we 

A-2 

1/2 
have  employed  the  relation  h  (x)  =  (jt/2x)  '   H  i(x)  between  spherical  and 

X'  X't"  2 

cylindrical  functions. 

The  other  choice  is  based  on  the  relation 

T^{-   cos  0)  =  (-l)^P^(cos  0),  (2.10) 

which  holds  for  integral  i,  and  leads  to 

♦{r,e)  =  -  iC  f{x.p,„)  P^.i  (-  cose)e""/==  ^^     .  (2,11) 


B.   Poisson's  Sum  Formula 

We  shall  now  consider  an  alternative  transformation  ,  based  on 
Poisson's  sum  formula  f 7) ,  which,  for  our  purposes,  may  be  written  as 
follows: 


°°  00        /-> 

Ilf(M)  =11  {-!)''[    f(x)  e 
£=0  m=-oo      -"o 


2imnX 


dX. 


(2.12) 


-7- 


Applying  this  to  (2.2),   we  get 

00 

00  r\ 


^     \{  COS  d)  exp 

A-2 


irtX(2m+  2") 


(2.13) 


XdA. 


Substituting  X  by   -X   in  the   integrals   for  m  =    -1  to   -  oo  and  making 
use   of  the  properties 


P  ^    \    (cos   Q)    =  P,     i(cos   9), 

-X-2  X-2 


„      (1)/    N  irtX  „   (1),    V      „     (2)/    N  -iitX^   (2)/    N 

H_^'    '(x)    =   e  H^'    '(x),    H_^'    '(x)    =e         H^'    '(x) 


(2.11|) 


(2.15) 


we  find  that  (2.13)  can  be  rewritten  as 


t(r,0)  =  C  V'^^' 
m=  0 


where 


fjr,9)    =  2(-l)'^£^  f(x,3. 


p)P^  i(cos  9)    exp 

X-2 


inX(2m+2) 


(2.16) 


XdA.    (2.17) 


This  result  is  equivalent  to  (2.7).   In  fact,  according  to  (2.1^4-) 

and  (2.15),  the  integrand  of  (2.7)  is  an  odd  function  of  X,  so  that  the 

lower  half  of  contour  C  may  be  replaced  by  its  reflection  about  the  origin, 

shown  in  broken  line  in  Fig.  1.   Accordingly,  contour  C  is  equivalent  to 

the  straight  line  D  located  above  the  real  axis.   On  D,  the  following  expansion 

is  valid: 

-iXK/2 


cos  nX 


=  2  e 


iXjt/2 


E  (-l)"e 
m=0 


m  2imjtX 


(2.18) 


-  8  - 

Substituting  this  result  in  (2.7),  we  are  led  to  (2.l6).   It  will  be  seen 
later  that  (2.l6)  is  a  more  convenient  form  for  several  purposes. 

C.   The  Incident  Wave 

We  shall  now  apply  Watson's  transformation  to  the  incident  wave  (2.1). 
The  resulting  expressions  will  play  an  important  role  later  on. 

A  derivation  similar  to  that  of  (2.7)  and  (2.11)  leads  to  the 
representations 

^ipcose  ^j^_^j  V2  ^_,^/^  y^  ^^^^^   G^(X,e)dA,  (2.19) 


ip  COS  ey^NV2  ^.i,/i,  ^       ^^^       ^^^^^  ^^  ^^^^^^ 


where 


2p    /  ^  J^   '^X''"    ^2' 


■iXjt/2 


G,(X,0)    =   X    e"       /      P,      i(cos   0)/cos   nX,  (2.21) 

X  A-    2 

G^(X,0)    =    -iXe^^'^/^   P,     i(-   cos0)/cos   jtX.  (2.22) 

d  X-  2 

The  next  step  is  to  deform  contour  C  into  a  contour  that  is  symmetrical 
about  the  origin.   For  this  purpose,  we  must  find  the  asymptotic  behavior  of 
the  integrand  at  infinity  in  the  X -plane. 

The  asymptotic  behavior  of  the  cylindrical  and  Legendre  functions 
is  discussed  in  Appendices  A  and  C,  respectively.   According  to  (A. 6), 


9  - 


and,  according  to  (C.8), 


,1/2 


G-^(X,0)  «  (2X/jt  sin  Q)  '         exp 


iXlf 


-  9   +  i  ^ 


( I  Xl  ->  00  in  I  )  , 


-i/p 
(2X/n  sin  0)  /   exp 


GgCx^e)  -«  (2X/n  sin  e)"'-/^  exp 

l/? 
«*  -(2X/jt  sin  9)  /   exp 


-ix(f -e)-if 


i^l-^  -if 


iM|-0J.if 


(|X|  ->«>  in  I  )  ,  (2.2l|) 


(  I  X I  ->  00  in  I  )  , 


X|  ^  00  in  I_)  ,  (2.25) 


where  I  and  I   denote  the  upper  and  the  lower  half  of  the  X-plane,  respectively. 

It  follows  from  these  results  that  contour  C  may  be  freely  deformed  in 
the  right  half -plane,  except  possibly  along  directions  approaching  that  of  the 
imaginary  axis.  To  study  the  behavior  along  such  directions,  let  us  introduce 
the  notations 


o  =  i  exp  (is)  =  exp 


1  I  2  +  ^ 


(2.26) 


T)  =  e  in|  2X/ep| 


(2.27) 


and  let  us  consider  the  behavior  of  the  integrand  when  X  ->  +  a  oo  and 
simultaneously  e  ->  0  in  such  a  way  that  t]  approaches  a  constant  value. 

The  curves  X  =  a|x|  ,  rj  ->  -n/2  and  X  =  -a|x|  ,  t]  ->  jt/2  for  large  |x| 
are  shown  in  broken  lines  in  Fig.  2  and  Fig,  3-   According  to  Appendix  A, 
the  function  J,  (p)  approaches  zero  to  the  right  of  these  curves,  and  so 
do  the  integrands  of  (2.19)  and  (2.20),  so  that  it  suffices  to  consider 


10  - 


their  behavior  to  the  left  of  the  curves.   In  this  region,  (2.23)  to 
(2.25)  imply 


I  (sin  e)^/^|j^(p)G^(X,0)|-  exp 


X|(t)+  0) 


for  X  ^  0  °o, 


and 


I  (sin  e)^/^\j^{p)G^{x,t 


exp 


exp 


exp 


-|X|  (t]  -  0  +rt) 


for  X  ^  -a  00,  (2.28) 


x|(t)  -  0) 


for  X  ->  a  00^ 


-\x\{ti  +  e-  n)      for  X  ^  -a  00.   (2.29) 


It  follows  from  these  results  that  the  path  of  integration  in  (2.I9) 
may  be  deformed  at  will,  provided  that  it  stays  asymptotically  to  the  right 
of  the  shaded  regions  in  Fig.  2.   A  symmetric  path  of  integration  can  be 
found,  as  shown  in  Fig.  2,  provided  that  the  corresponding  t]   satisfies 

0-:t<ri<-9,  (2.30) 

which  is  only  possible  for  0  <  :rt/2. 

Similarly,  as  shown  in  Fig.  3^  a  symmetric  path  of  integration  can 
be  found  in  (2.20),   provided  that  the  corresponding  tj  satisfies 
jr  -  0  <  T)  <  0, 


(2.31) 


which  is  only  possible  for  0  >  n/2. 

We  shall  therefore  take  (2.I9)  with  the  path  of  Fig.  2  for  0  <  n/2, 
and  (2.20)  with  the  path  of  Fig.  3  for  0  >  n/2.   Making  the  substitution 

X  -»  -X  in  the  integrals  from  -a  00  to  0  and  employing  (2.1^+)  as  well  as  the 
identities 


imX 


rj-^Q.w 


'ReX 


V^f 


TT 
2 


Fig.    2 


7^^^     (Too 


ReX 


Tj-^TT-Q 


Fig.    3 
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J^(p)  -  e^'^^J_^(p)  =  -  ie^'^^  sin  nX   h[^^(p),  (2.32) 

J^(p)  -  e'^'^^J.^Cp)  =   ie"""--  sin  nX  i^^kp)  >  (2.33) 

2 
we  finally  get  the  integral  representations 


O  00 

ip  COS  e^    J-2!-Y^^   e^«/^  ^     li^hp)    R  i    ( cos  9  )  e^W2tan(  nX)Xd> 
\  2d   /  .J^      X  X-2 


\    y  0  (2.34) 

(e  <  I), 

I  r-,  0  00 

^ipcos  e  J^y       e""/^  C     H[2^(p)P^_i(-cos  e)e"^^'^/2tan(KX)Xdl 

rt    (2.35) 

the  corresponding  paths  of  integration  being  those  shown  in  Fig.  2  and 
in  Fig.  3 J  respectively. 


12 


III.   The  Poles  of  the  S -function 

In  order  to  deform  the  path  of  integration  in  (2.7)  and  (2.11)  away 
from  the  real  a:xis,  we  need  information  about  the  singularities  of  the 
integrand  in  the  X-plane.   The  integrand  is  a  meromorphic  function  of  X, 
and  its  poles  are  the  poles  of  the  S-function 

S(X,3)  =  -H[2)(p)/H[l)(f3),  (3.1) 

which  are  the  roots  X  ( p)  of  the  equation 

h[^^|3)  =  0.  (3.2) 

They  might  be  called  the  Regge  poles  for  the  hard  sphere  problem,  although 
they  do  not  show  the  typical  Regge  behavior  characteristic  of  Yukawa-type 
potentials. 

The  roots  of  (3-2)  have  been  discussed  by  several  authors  (2,^ ,9 ,10) • 
We  are  only  interested  in  their  behavior  for  large  values  of  p. 

It  follows  from  (2.15)  that  the  roots  are  symmetrically  distributed 
with  respect  to  the  origin,  so  that  it  suffices  to  consider  them  in  the 
right  half -plane.   In  this  region,  there  exists  an  infinite  number  of  roots, 
all  located  in  the  first  quadrant,  close  to  the  curve  h  defined  by  (Cf.  Appendix 

Re[(X^  -  3^)^/2  -  X  cosh"^  (x/p)]  =  0.  (3-3) 

This  curve  is  shown  in  broken  line  in  Fig.  ^.   It  cuts  the  real  X  axis  at 
X  =  3,  at  an  angle  of  Tt/3,  and  tends  to  become  parallel  to  the  imaginary 
axis  as  |x|  -*  oa.      All  the  roots  are  simple  (£) . 

The  roots  of  greatest  physical  importance  are  those  closest  to  the 


ImX 


^ReX 


Fig.   h 
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real  axis,   which    are    located  in  the  neighborhood  of  X  =  p.   In  this 
region,  we  can  use  the  expansion  (A. 17): 


K 


1) 


X  -X  e^"/3  (x/2)l/3 


2  e 


■^"/3(2/x)V3Ai(-x)  .^(x-1).  ^3.^) 


Since  the  zeros  of  the  Airy  fimction  are  all  located  on  the  negative 
real  axis,  we  get  (k) 


X^(p)  =  3  +  (3/2)^/3  ^^^U/3  ^   ^p-1/3)^ 

where  -x  is  the  nth  zero  of  the  Airy  function.   A  table  of  the  first  five 
zeros  is  given  in  Appendix  D.   According  to  (D.7),  we  have,  for  large  n, 

^/3  U/3  a/3 


(3.5) 


^Jp)'*  P  +  2 


3n(n  -  y-  ) 


P 


(n  »  1) 


(3.6) 


Let  us  finally  consider  the  asymptotic  behavior  of  the  roots  for 
|x|  -^00.   In  this  region,  the  behavior  of  H;   in  the  neighborhood  of  h  is 
given  by  (A.7) : 


h[1)(p).  2(2/.)V2(,2_32^-VU^i.A^i^ 


(x2_p2)l/2_,^^ 


so  that   the  roots   are  given  by 


X   £n 
n 


X   .   (x2-p2)V2 
n         n 


,,2    ^2,1/2        .,  1   > 

-   (X„-3    )    '       =    i(n    -  ^  )    jT. 


,       ,2   „2>l/2\ 

X   +{X    -P    )    M  rt 

p        y   '  ^ 

(3.7) 


(3.8) 


Let 


P^   exp 


'2    "   ^n^ 


(3.9) 


-  ll+  - 


Then,  (3-8)  gives 


e  in(2p  /ep)«  it/2, 


where 


This  equation  may  be  solved  by  iteration.   We  find 

-1 


p^  =  (n  -  ^)  n  hn   2jt(n-  ^)/e3 


1  +  & 


^n^n  n 
in  n 


(3.10) 
(3.11) 


(3.12) 


so  that  (10) 


X  (3) 


-lit 


^"("  -  it)     exp  /     ^ 

^-^^-"1      I  2  in[2«(n-^)/e3] 


(3.13) 


?n[2rt(n  -i-)/e3] 


1  +  <yi 


in  in  n 
,  in  n 


Thus,  both  Re  X  and  Im  X   approach  infinity  with  n,  but  Re  X   does  so  more 
'         n        n  ^  ^  ^ 


(|^n|»  P  )• 
-It 


slowly  by  an  inverse  logarithmic  factor.  Note  also  that  X   -  X  ^=  ©'[(in  n)   J, 

so  that  the  roots  cluster  closer  together  as  n  ->  oo. 

The  residues  of  the  S-function  at  its  poles  will  also  be  required. 
According  to  (3.I), 


r  ( 3)  =  residue  S(X,3) 


X=X 


=  -h(='(3,/^h(^)(3, 


n     n 


X=X 


(3.1i+) 


In  the  neighborhood  of  X  =  3  we  have,  by  (A.I7), 


H. 


(2) 


X  -X   e^"/3(x/2)l/3 


=  2  e^^/3(2/x)^/3  Ai(-x  e^WS)  ,  ^^-1),   ^^_^^^ 
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It   follows   from  (3-1^),    (3-15),    (3-^)    and  (D.2)    that 


r  (3); 
n 


2n  ^2;  [Ai'(-x„)]2 


:3.i6) 


at  the  poles  (3- 5)-   Asymptotically,  for  n  »  1,  we  may  employ  (D.8),  which 
leads  to 


^  1  -in/6  1/3 
r  ~  —  e   '  3  '  -" 
n   2 


3n(n  -  I  ) 


-1/3 


(n  »  1). 


(3.17) 


Finally,  let  us  consider  the  residues  at  the  poles  (3.13)-   In  this 


(2) 
region,  }r    is  given  by  (  cf.  Appendix  A) 

h[2)  (3)-  i(ff  V^-3^)-^/^  expi-(X^-32)V2  ,  ,  ,, 


(3.18) 

/,2  ,2, 1/2-1' 
X  +  (X  -3  )  ^ 


so  that,  according  to  (3.8), 

n 
Similarly,  according  to  (3-7)  and  (3.8), 


(3.19) 


a4^^3)/ax 


x=x 


(-l)-l2e-/\2/.)V2(,2.p2^-l/U 
2    2,  1/2 


in 


X  +(X   -  3 
n   n 


It  follows  from  these  results  that 


(3.20) 


rj^(3)=  |j-2n[2«(n-  ^  )/e3]    \       exp 


lit 


2in[2jt(n   -  ^)/epJ      , 


(\\\»n, 


(3.21) 


so  that  r  ->  0  like  (in  n)    for  n  ->  00. 
n         ^    ' 
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IV.   The  Residue  Series 

Let  us  consider  the  integral  th^  appears  in  (2.17).   We  already 
know  that  the  integrand  has  an  infinite  niimber  of  simple  poles  in  the  upper 
half -plane.   Let  us  now  inquire  under  what  conditions  the  integral  is  reducible 
to  a  series  of  residues  taken  at  these  poles.   For  this  purpose,  we  must  find 
a  sequence  of  paths  C  passing  between  the  poles  and  such  that 


lim  \   f(X,3,p)P^  _i(cos  9)    exp 


n  -^  00  "^C 


inX(2m  +  i) 


XdX   =  0. 


(h.l) 


Let  us  consider  the  behavior  of  the  integrand  as  |x|  ^oo  in  I  .   It 
is  shown  in  Appendix  B  that 


^(X,P,p)«  ^ 


f) 


X 


(|X|   ->c»), 


(1^.2) 


:(1) 


where  g  is  given  by  (2.9). 

According  to  Appendix  A,  the  behavior  of  H^"^''(3)  differs  on  the 
right  and  left  of  the  curve  h  (cf.  (3-3)  )  where  the  poles  are  located.  In 
regions  2  and  3  of  Fig.  h^   we  have 


whereas,  in  region  1, 


H<^'0) 


2  \l/2  ^^a 


nX 


2X/   ' 


1/2  /2X\^ 
eS 


(^.3) 


(U.U) 


In  the  neighborhood  of  h  ,  according  to  (A. 7),  we  must  take  the  sum  of  both 
estimates: 


"!"<«  »^,x 


sinh 


X  t 


"(i) 


IT 


(U.5) 
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Let  us  exclude,  for  the  moment,  the  directions  9=0  and  0  =  n. 
Then,  according  to  (C.8), 

V     j^(cos9)  «  (2jTXsin0)~"'"/^  exp( -iX9  +  irf)   (Ixl^ooinl). 

It  follows  from  these  results  that,  for  |x|  -^  oo,  the  integrand  of 
(U.  1)  behaves  like 


(U.6) 


exp 


exp 


exp 


X^n{  ^)+  iX  f|  -  0]+  2imrtX 


X  tn  - 


2Xr 


2 
\eka  / 


-V  ix(| 


9U-   2imjtX 


-X  enfS^L' 


iX(  I  -  9)+  2imjrX 


in  region  3> 


in  region  2, 


in  region  1, 


(U.7) 
(U.8) 

(h.9) 


except  in  the  neighborhood  of  curve  h  ,  where  it  behaves  like 

^n(r/a)  +  i^  |  "  ^+   2imjtX 


K(X)  = 


exp 


sinh 


C+.io) 


X  en(2X/ep)  +15- 

It  is  readily  seen  that  (^.7)  approaches  zero  at  least  exponentially 
for  all  ©,  0  <  0  <  n;    the  same  is  true  for  (h.Q)    and  {h.g)    if  m  >  1.   However, 
for  m  =  0,  (^.8)  and  (^.9)  go  to  zero  everywhere  if  and  only  if  0  <  it/2- 

Finally,  near  h  ,  we  have  to  avoid  the  poles,  which  are  the  zeros 

of  the  denominator  of  (^.10)  (of.  (3.8)  ).   For  this  purpose,  we  shall  choose 

r  as  a  half -circle  of  radius  R  passing  half-way  between  consecutive  poles,  so 
^n  n  ^ 

that,  for  large  n,  we  have,  according  to  (3.II), 

R^in(2R^/e3)  +  f  =  nn  +  |  .  ('^•ID 
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Then,  in  the  neighborhood  of  h  , 


X  =  R  exp 

n 


i  (|-  e) 


:i^.i2) 


and,  since 


2     2         2      '-'   1/2 
sinh  (a  +  ib)  I  =  ( sinh  a  cos  b  +  cosh  a  sin"lD)  '  >  cosh  a  |  sin  b|  , 


it  follows  from  (^+.10)  that 


(U.I3) 


K(X)|  < 


exp 


-R  (2m7r  +  ^ 
n       2 


^n  -  ) 

a  , 


< 


cosh  (R  [e  in(2R  /ep)  -  ^] 


<  exp[-R^(2m«  +  |  -  6  -  e  tn  |  )J  . 


(i^.lil) 


Since  e  -> 0  along  h  and  its  neighborhood  ,  we  see  that  K(\)  approaches  zero 
exponentially  for  0<9<n,  m>l  and  for  0  <  9  <  7t/2  if  m  =  0. 

Since  P   i(l)  =  1,  it  is  readily  verified  that  all  the  above  results 
remain  true  if  0  =  0.   However,  near  0  =  n,  we  can  no  longer  employ  (2.17), 
since  P  i(cos  9)  has  a  logaritlimic  singularity  at  this  point. 

In  conclusion,  we  see  that  (^.1)  is  valid  for  all  m  >  1,  0  <  9  <  n, 
but  it  is  only  valid  for  0  <  9  <  Jt/2  if  m  =  0.  It  then  follows  from  (2.1?) 
that 

2p" 


t  (r  6)  =  (-l)"'2ne""/^  ( 
m  \ 


m„   ijt/^  (    IT 


•C  ^n^n  exp[iTTX^(2m  +  \)\    h[^^(p)P^   _   1    (cos   9) 
n=l  -      ? 


provided  that 


(^.15) 


m  =   0,    0  <  9  <  ^     or  m  >  1,      0  <  9  <  jt, 


(i^.l6) 


19 


where  X   ( p)  are  the  poles  of  S(X,3)  and  r  ( p)  is  defined  by  (3.1^). 

One  can  verify  directly,  with  the  help  of  (3. 13)  and  (3. 21),  that, 
for  n  -^  00,  the  nth-order  term  of (^.15)  behaves,  in  absolute  value,  like 


1  njT 

^"^    "  2tn(2nn/ep) 


o     ^       a  Jtin(r/a) 


so  that  the  conditions  for  the  convergence  of  the  residue  series  are  exactly 
those  stated  in  (^.l6). 

A  rigorous  discussion  of  Watson's  transformation  and  the  convergence 
of  the  residue  series  in  the  shadow  region  was  apparently  first  given  by 
Pf lumm  ( 11) . 
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V.  The  Shadow  Region 

We  shall  begin  the  discussion  of  the  solution  with  the  simplest 
case,  namely,  the  behavior  in  the  shadow  region. 

According  to  geometrical  optics,  the  shadow  region  is  the  whole 
cylinder 

0  <e  <  e     <  n/2,  (5.1) 

where 

e     =  sin"^(a/r)  (5-2) 

o  ^  /   ' 

is  the  shadow  boundary  angle.   Actually,  as  will  be  seen  later,  the  shadow  of 
the  sphere  does  not  extend  beyond  distances  ~  ^a,    and  the  transition  to  the 
illuminated  region  already  starts  at  much  smaller  distances,  of  the  order  of 
pl/3a. 

Since  9  <   n/2  in  the  shadow  region. (^. 15)  is  valid  for  all  m  >  0,  so 
that  (2.16)  becomes 


Hr,e)    =2ne"/V^f^'    •      £       (-1)" 
V^P/  m  =   0 


00  ,         , 

•J"'     X   r      exp    [ijtX    (2m  +  |)]    Hi"^^(p)P,       1    (cos    9 
n  =  1  '■^  n 


.     /I,  ,         •>l/2    00  exp( -ijtX  /2)         ,^. 

,eWn^]'      YZxr  — ^ H^^P)P,      i(cos0),        (5.3) 

V  2p/  '— »  ^    n  n  cos   txX  X„      "^     X„-2  \-^   ->/ 

•^'n=l  n  "  n^ 


where  we  have  employed  (2.l8).   The  last  form  of  this  result  could  also  have 
been  obtained  directly  from  (2.7). 
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Eq.  (5.3)  is  formally  similar  to  an  "eigenfunction"  expansion  in  terms 

of  the  functions  Hj   (p),  which  satisfy  the  boundary  condition  on  the  surface 

n 
of  the  sphere  (cf.(3-2))  and  the  radiation  condition  at  infinity.   The  latter 

implies  a  non-self adjoint  problem,  explaining  why  the  "eigenvalues"  are  complex. 
It  was  remarked  by  Sommerfeld  (12)  that  the  "eigenf unctions"  are  even  orthogonal 
in  a  certain  sense,  and  he  proposed  to  derive  (5-3)  on  this  basis.   However, 
as  was  shown  by  Pflurran  (ll),  the  set  of  "eigenfunctions"  is  by  no  means  complete, 
and  no  general  characterization  of  the  class  of  functions  for  which  the  expan- 
sion is  applicable  has  so  far  been  given. 

In  practice,  even  though  (5-3)  converges  for  all  Q<   n/2,  its  usefulness 
is  restricted  to  the  domain  where  its  terms  are  rapidly  decreasing  from  the 
beginning,  so  that  only  the  first  few  terms  have  to  be  considered.   This 
happens  only  within  the  shadow  region,  as  will  now  be  seen. 

Let  us  consider  the  first  few  terms  of  (5-3)j  corresponding  to  poles  of 

the  type  (3-5)j  located  near  X  =  p.  We  shall  restrict  ourselves  to  points 

within  the  shadow  region,  not  too  close  to  the  surface  of  the  sphere,  so 

that 

r  -  a  »  p"^/3  ^_  (5.I1) 

Under  these  conditions,  we  have  kr  -  \\    |  »|  >^j^|    >    so  that  we  may  employ 
the  expansion  (A.  16)  for  h;  ^(p).   Furthermore,  ass\iming  that 

e  »  p"^,   ""  (5.5) 

we  have  Ix  |0  »  1,  so  that  we  may  employ  expansion  (C.8)  for  P   ^(cos  9). 
'  >^'  n  2 

Substituting  these  expansions  in  (5-3),  and  taking  into  account  (3-5) 
and  (3.16) ,  we  get 
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e 


-lit 


/6  /fi\V3    /  2         \l/k 


Uy      I   2      2   i         ~~     :    ,,a/2 


^^^''^%,  a/2       I2;      l^ ^i  ;:    '■    ,,1/2 

(2jt)  \r      -   a      /  (kr   sm  9)    ' 


m     ^— 7  1 


'^n    ^"^^"^     ^    FTT; ^72    [^^^i  Vm  +   i  f)    +   exp   (i  X^5    -i  f  )]  , 

m  =  0  n       [_Ai '  (  -x^)  ] 

(5.6) 


where 


7      =9      -  0   +   2imT,  (5.7) 

'mo  ' 

6      =   e     +  0  +   2nw  (5.8) 

mo  \  -^      / 

and  9   is   given  by  (5-2). 

In  writing   (5-6),   we  have   already  assumed  that   only  the   first   few  terms 

of  the   series   give  a  significant   contribution,    since  the   approximations   employed 

correspond  to  the  poles    (3- 5)-      Thus, 

1/3 
|exp(iX^7„)|    =   exp    [    -  ^x^  [|j  (0^   -   0)]  (5-9) 

must  be  a  rapidly  decreasing  function  of  n.   This  will  be  true  provided  that 

7^  =  ^O  -  9  »  f^'^,  (5.10) 

so  that  (5-6)  is  valid  in  this  region. 

The  physical  interpretation  of  (5-6)  is  well  known  (13.)-   At  short 
wavelengths,  one  may  employ  the  concept  of  propagation  along  rays.   The 
incident  rays  that  are  tangent  to  the  sphere  at  T-j^  and  Tp  (Fig.  5)  excite  a 
series  of  surface  waves  emanating  from  thes'e  points.   These  waves  travel  along 
the  surface  with  phase  velocity  slightly  smaller  than  that  in  free  space,  due 
to  the  delay  in  overcoming  the  curvature  of  the  sphere.   As  they  travel  along 


(r2.a2)'/^ 


P(r,^) 


(r2-a2)'/2 


Fig.    5 


-  23  - 

the  surface,  they  shed  radiation  along  tangentis    '  ections,  leading  to  the 
angular  damping  factor  (5-9) • 

A  point  P  within  the  geometrical  shadow  is  reached  by  two  rays  emanating 
from  the  points  T'  and  T'  where  the  tangents  to  the  sphere  from  P  meet  the 
surface  (Fig.  5)-   The  corresponding  angles  travelled  along  the  surface  are, 
according  to  Fig.  5^ 

7^  =  I  -  cos"\a/r)  -   9  =   9^   -   Q, 

5  =  ^  -  cos"^(a/r)+  9=9     +9, 
o   2         \  /  /       o    ^ 

in  agreement  with  (5.7)  and  (5.8)  for  m  =  0.   The  paths  T  T'P  and  T  T'p  are 
called  "diffracted  rays"  in  Keller's  geometrical  theory  of  diffraction  (13)  . 

The  terms  with  m  >  1  in  (5-6)  correspond  to  rays  which  have  encircled 
the  sphere  m  times  before  leaving  the  surface,  so  that  the  corresponding 
angular  paths  are  increased  by  2mjT.   This  interpretation  is  corroborated  by 
the  solution  to  the  problem  of  diffraction  of  a  pulse  (lif),  where  one  can 
follow  the  diffracted  wave  front  around  the  sphere.   In  each  encirclement,  the 
rays  go  through  the  points  T,  and  Tp,  where  all  diffracted  rays  meet.   As  is 
well  known,  passage  through  a  focal  point  leads  to  a  phase  decrease  by  n/2. 
This  explains  the  factor  [exp  (-in/2)]    =  (-1)   in  (5.6). 

The  phase  factor 

exp  {ik(r^-a^)^/2  +  ik[l  +  x^(^)"^/3j  a  y^] 

also  corresponds  to  the  described  optical  path  (and  similarly  for  the  term  in 

5  ).   Note  the  decrease  of  the  phase  velocity  along  the  surface. 

m  ^ 

2  2  1/*+  1/2 

The  denominator  (r  -a  )    (r  sin  9)  '   can  also  be  obtained  from  the 

law  of  conservation  of  the  intensity  along  a  pencil  of  rays  (13)-   This 

'denominator  would  vanish  at  r  =  a,  which  is  a  caustic  of  diffracted  rays,  but. 
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according  to  {^.^) ,    (5.6)  is  not  valid  there  (the  actual  value,  of  course,  is 
t{a.,6)    =0).   It  would  also  vanish  at  0  =  0,  but  this  is  excluded  by  (5- 5). 

The  direction  9  =  0  is  a  focal  line  of  diffracted  rays.   In  fact,  an 
observation  point  P  on  the  axis  is  reached  by  a  whole  cone  of  diffracted 
rays  tangent  to  the  sphere,  instead  of  by  two  rays  only.   This  focussing  effect, 
which  is  responsible  for  the  well-known  Poisson  spot  (cf.  Section  VIII),  leads 
to  an  enhancement  of  intensity  near  the  axis. 

In  fact,  for  9  <  p"  ,  we  must  employ  (C.9)  instead  of  (C.8),  so  that 
(5.6)  is  replaced  by 

sl/3.  ,1/2/     2   Nl/l| 


t(r,e) 


00  expFiX    (2mjT+9  1 

m=0  n         [Ai'(-x  )] 


The  surface  waves  associated  with  the  poles  X  are  also  known  as 

n 

"creeping  modes"  (5).   It  should  be  emphasized,  however,  that  the  above 
physical  interpretation  applies  only  to  the  lowest -order  modes.   As  soon  as 
the  damping  within  a  single  wavelength  becomes  appreciable,  the  above  concepts 
lose  their  validity.   In  practice,  of  course,  only  the  lowest-order  modes  give 
a  significant  contribution. 
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VI.   The  Nearby  Lit  Region  and  the  WKB  Approximation 

We  shall  now  consider  the  behavior  of  the  solution  in  the  lit  region, 
not  too  close  to  the  surface,  so  that  (5-^)  is  assumed  valid,  but  still  not 
in  the  far  field  region,  which  will  be  treated  in  Section  IX.   We  shall  also 
stay  away  from  the  neighborhood  of  the  shadow  boundary,  which  will  be  dis- 
cussed in  Sections  VII  and  VIII. 
A.   The  Forward  Half-Space 

For  0  >  0  ,  we  have  7  <  0  in  (5.6).   According  to  (5-9),  this  implies 
that  the  residue  series  containing  7  starts  out  with  exponentially  increasing 
terms.   It  does  not  follow  that  the  residue  series  representation  becomes  in- 
correct, for,  as  we  have  seen  in  Section  IV,  it  still  converges  for  9  <   n/2. 
What  happens,  however,  is  that  it  becomes  useless  for  all  practical  purposes. 
The  physical  reason  for  this  behavior  is  that  the  wave  function  is  no  longer 
exponentially  damped,  but  contains  additional  contributions  corresponding  to 
the  incident  and  reflected  waves  in  geometrical  optics. 

A  modification  of  Watson's  transformation  to  take  into  account  these 
contributions  in  the  lit  region  was  first  proposed  by  Fock  (3.)  •   As  reformulated 
by  Franz  (J^>5)7  the  basic  idea  is  to  substitute  in  (2.7),  taken  over  the  alterna- 
tive contour  D  in  Fig.  1,  the  identity  (C.5) 


P  j|^(  cos  9)  -   -  le 
X-2 


P   i(-cos  9)    +   2i  cos  itXQ;  i  (cos  9) 
X -2  X -2 


(6.1) 


The  integral  of  the  first  term  on  the  right-hand  side  of  this  expression  is 
then  reduced  to  a  residue  series,  which  gives,  as  we  shall  see  later,  the 
continuation  of  (5-3)  into  the  lit  region.   The  second  term  cancels  the 
denominator  cos  jtX  in  (2.7),  so  that  the  integrand  no  longer  has  poles  at 
the  positive  half -integers,  and  the  integral  can  be  evaluated  by  the  saddle- 
point  method,  yielding  the  contributions  from  the  incident  and  reflected  waves. 
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We  shall  see,  however,  that  the  proposed  contour  of  integration  can  only  be 
employed  in  the  forward  half -space  (9  <  n/2).      In  this  region,  we  shall 
derive  essentially  the  same  results  by  a  different  method,  which  has  the 
advantage  of  greater  simplicity,  as  well  as  of  showing  more  clearly  the 
connection  between  the  lit  and  shadow  regions.   This  method  does  not  require 
the  re-evaluation  of  the  whole  residue  series,  but  only  of  that  part  of  the 
term  m  =  0  in  (5.6)'  that  "goes  wrong"  for  6  >  9   . 

In  fact,  for  m  >  1,  not  only  is  the  transformation  that  led  to  the 
residue  series  allowed  for  0  <9  <  n   (cf.  (^.l6)),  but  also  the  corresponding 
residue  series  in  (5-6)  are  all  rapidly  convergent,  so  that  we  need  only  be 
concerned  with  the  term  m  =  0.   The  corresponding  term  in  (2.17)  is 

00 

t  (r,0)  =  2  r  f(X,3,p)P,  i(cos0)e^^"/^  X6X.  (6-2) 

Let  us  make  the  splitting  (cf.  (C.2)) 

P,  i(cos  0)  =  q!^H    (cos  9)    +  J^l  (cos  9).  (^-3) 

A'~2  ^  ~2  -^"2 

Then,  according  to  (C.?),  the  term  in  y      in  (5.6)  arises  from  Q._i  ,   while 

(2) 
that  in  5  arises  from  Qr    \- 

The  residue  series  containing  &  remains  rapidly  convergent  in  the  lit 
region,  provided  that  (cf. (5-8)) 

3V3  (0  +  0)  »  1,  (S.k) 

o  ' 

which  will  be  ass\imed  throughout  this  Section.   The  vicinity  of  the  forward 
direction,  where  this  assumption  is  not  fulfilled,  will  be  considered  in 
Section  VIII. 
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Thus,   we  finally  have  to  consider  only  the  expression 
00  +  ie 


i^^(r,e)    =  2^    _  f(X,p,p)Q[H  (cos  9)  e^^''^^^,  (6.5) 


where  the  path  of  integration  has  to  be  taken  slightly  above  the  real  axis 
in  order  to  avoid  the  poles  of  q}  ^  at  the  negative  half-integers.   This  is 
done  before  the  splitting  (6.3).   We  shall  see  that  the  geometrical-optics 
approximation  to  the  solution  is  entirely  contained  in  this  integral. 

The  behavior  of  the  integrand  of  (6.5)  as  |x|  ^  co  in  I   is  similar 
to  that  of  (^.1)  with  m  =  0.   In  fact,  according  to  (^.6)  and  (C.7),  the 
behavior  of  P  j^(cos  9)    in  I   is  the  same  as  that  of  Q;  j^  (cos  9).      The 

A-2  +  ^~2 

discussion  given  in  Section  IV  implies  that  the  path  of  integration  in  (6.5) 
may  be  deformed  at  will  in  1  ,    for  9  <   jt/2. 

Following  Franz's  method  (_5)  ,  we  shall  deform  it  into  the  path  P  shown 
in  Fig.  6,    going  around  the  poles  of  S(X,3),  and  beginning  and  ending  at 

infinity  in  the  region  between  the  curves  t)  ->  3n/2  andr]-*  -jt/2,  where  r)  is 

(2) 
defined  by  (2.27).   According  to  Appendix  A,  H)  '(p)  -^  0  as  |x|  ->  00  in  this 

region.   It  follows  that  we  may  split  (6.5)  into  two  integrals,  corresponding 

to  the  two  terms  in  the  second  member  of  (2.8),  since  each  of  the  integrands 

will  separately  go  to  zero  as  |x|  ->  00  in  this  region.   The  integral  containing 

(2) 
YL      (p)  identically  vanishes,  since  the  contour  may  be  closed  at  infinity  and 

the  integrand  has  no  singularities  within  i't.   Finally,  we  get 

*i"'^.^)  =  -(§?/'  e-'"/^  J^  ^    4"(P)  i%   (cos  e,e'^"/=,^ 


1^(3) 


(e  <  h/2). 
(6.6) 


ImX 


(7200  77—-^ 


rReX 


\=psin9 


—  TT 


Fig.    6 
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The  path  T   crosses  the  real  SLxis  twice:   first  between  0  and  3 
and  then  "between  3  and  p.   We  shall  see  that  there  is  one  saddle  point  of 
the  integrand  in  each  of  these  intervals,  and  P  will  be  taken  through  both 
saddle  points.   In  the  lower  half -plane,  to  the  left  of  the  curve  t\  -*  n 
(Fig.  6),  the  integrand  increases  exponentially  for  |x|  -» oo.   However,  the 

steepest  descent  paths  corresponding  to  the  two  saddle  points  can  be  joined 

(2) 
by  an  arc  going  through  the  neighborhood  of  the  first  zero  of  H/   ( 3) ;  where 

the  integrand  is  small  (5)  • 

To  the  right  of  the  curves  h  and  h  in  Fig.  6,  we  have,  according 

to  Appendix  A, 

h[2)(p)/h[i)(3)  «  -  1.  (6.7) 

Furthermore,  we  can  employ  the  expansion  (A.l6)  for  K   (p)  and  the  expansion 
(C.7)  for  q}  {    (cos  9).   Making  the  change  of  variable 

X  =p  sin  w,  (6.8) 

we  find  that  the  portion  of  the  contour  T   to  the  right  of  h  and  h  contributes 

1/2 
*1=   (  g/sin  9)       ""'"'^^  yA(w,p,0)  exp  [ipa(w,0)]  dw,   (6.9) 

where  the  path  of  integration  crosses  the  real  w  axis  between  w  =  0  and 
w  =  n/2,  and 

a(w,0)  =  cos  w  +  (w-0)  sin  w,  (6,10) 

A^  a\         r    ■  n1/2    C  i      /cot   9  3+2   siQ^  w\ 

A(w,p,0)    =   (sin  w  cos  w)   /        1  +     I^  (   2   sin  w     "     , SF 

{  ^    \  6  cos     w  / 

+  »(p"^)]-    (6.11) 
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In  this  approximation,  (6,9)  is  independent  of  the  radius  of  the  sphere, 
so  that  ^,  should  correspond  to  the  incident  wave.   This  is  indeed  so,  as 
will  now  be  seen. 

The  exponent  (6. 10)  gives  rise  to  a  saddle  point  at 


corresponding  to  \  =  p  sin  Q ,    as  shown  in  Fig.  6.   The  corresponding  steepest 
descent  path  crosses  the  real  axis  at  an  angle  of  rr/^. 

Taking  into  account  not  only  the  main  term,  but  also  the  first  correction 
term  in  the  saddle-point  method,  the  behavior  of  (6.9)  for  large  values  of 
p  is  found  to  be  given  by 


^ 


Ae 


ipa 


^   (a"  sin  0) 


1/2 


1  + 


2a"p 


A"  A'  g/1   l_/a:lr   1  ^    , 
A  "a  a"  ^  12  l^a"  J      "  IT  a"  J 

+  er(p-2)j,     (6.12) 


where  A,  a  and  all  their  derivatives  (denoted  by  A',  a",  ...)  are  to  be 
taken  at  the  saddle  point . 

Substituting  A  and  a  by  (6.10)  and  (6.11),  it  is  found  that  the 
expression  within  square  brackets  identically  vanishes,  so  that 

(p  »  1). 


^  ^iPcos  0^  ^^p-2^ 


(6.13) 


Thus,  the  contribution  from  the  right-hand  saddle  point  is  essentially 

identical  to  the  incident  wave. 

To  the  left  of  the  curves  h  and  h  (Fig.  6)  we  have,  according 
to  (A.  16), 
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.(2),„^/Jl), 


h;^^^P)/H^-'^(P)  =  exp<'-2i  (3  -n^/  -  X  cos-^(x/3)-  f 


-1, 


1  + 


hTW^ 


1.?  ^ 


3    n2      .2 

■^  p  -  X 


i^.lh) 


and,  again  making  the  change  of  variable  (6.8),  we  find  that  the  contribution 
from  the  portion  of  the  contour  P  to  the  left  of  h  and  h  is 


^  = 
r 


2n   sin  9 


V2  u/k 


j  B(w,p,0, 


7  J  exp 


ip&(w,0,7] 


dw, 


(6.15) 


where 


7  =  a/r. 


(6.16) 


S(w,0,7)  =  cos  w  -  2(7   -  sin  w)  ' 


+  sm  w 


w  -  0+  2  cos 


-1  f  sin  w 


(6.17; 


B(w,p,0,7)  =  (sin  w  cos  w  ) 


V2  li . 


I? 


2      2 
cot  0      (2sin  w  +  37  ) 

2  sin  w  "^  77"2     .  2.3/2 
L-  3(7   -  sm  w)  ' 


or 


(3+2  sin  w) 


D  cos  w 

The  saddle  point  is  determined  by  the  condition 

.  -        /e  -  w 
sm  w  =  7  cos  I  — — 


+   &ip' 


(6.18) 


p  =  r  sm  w  =  a  cos 


w 


(6.19) 


The  corresponding  value  of  X  is  X  =  p  sin  w  =  kp  (cf.  Fig.  6).   The  steepest 
descent  path  crosses  the  real  axis  at  an  angle  of  -  n/h. 


P(r,^) 


w7=^-2£ 


f-C 


?-v^ 

1            / 

p/           \ 

0 

Fig.    7 
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The  physical  interpretation  of  (6.19)  is  shovm  in  Fig.  7.   According 
to  (2.5),  p  =  x/k  may  be  interpreted  as  the  impact  parameter  associated 
with  an  incident  ray.   As  shown  in  Fig.  7?  this  is  precisely  the  incident 
ray  AB  that  reaches  the  observation  point  P  after  being  reflected  at  the 
surface  according  to  the  laws  of  geometrical  optics.   The  angle 

^=|(0-w)  (6.20) 

is  the  complement  of  the  angle  of  incidence. 

A  saddle-point  evaluation  of  (6.I5),  including  the  first  correction 
term,  yields,  similarly  to  (6.12), 


t     = Be 


ip6         (       .     7-Dii  -Q'j^'"    •-      /c-,,,\    -,   .,IV 


r     /  I  ciii  __•..  ON  1/2 


B"  B'&"'    ^  f^"'l   1  2_ 
B  ^B  T&^l"^  12  V&"  ,/    ^  1&"| 

-20   (6.21) 


(l&"lsin  9)"/"   I    2p|5"| 

+  &(p"^)]  . 

where  B,S  and  their  derivatives  are  to  be  evaluated  at  the  saddle  point. 

It  is  convenient  to  express  the  result  in  terms  of  the  parameters 
^  aind 

s  =  r  cos  w  -  —  sin  ^,  (6.22) 

which  measures  the  distance  (taken  along  the  ray)  from  the  observation  point 
to  the  caustic  of  the  reflected  rays.   Substituting  (6. I7),  (6. 18)  and  (6. 19) 
in  (6.21),  it  is  found  that 
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r 


2  .   o;. 
a  sm  2^ 


-  *+s(s  sin  2  ^  +  a  cos^' 


1/2 


exp 


ik(s  -  ^  a  sin  ^) 


1+ 


2p 


sin^^ 


~T-        2     2~~ 
2  sin  ^  cos  ^ 


s    5  V^^"^^ 


f-^(:f;r,-53ind(f 


15    2,  /a 


2  sin  ^  cos  ^    (s  sin  2^+a  cos  ^) 


hC(k 


_2l  (o-23) 


This  asymptotic  expansion  in  inverse  powers  of  k  corresponds  to  the 
well-known  WKB  approximation,  which  has  been  investigated  by  Luneburg  (l5) 
and  Kline  (l6)  in  connection  with  Maxwell's  equations. 

The  first  term  of  (6.23)  represents  the  reflected  wave  according  to 
geometrical  optics  (first -order  WKB  approximation).  The  amplitude  of  this 
term  takes  into  account  the  divergence  of  the  rays  after  reflection  at  the 
surface. 

The  remainder  of  (6.23)  represents  the  correction  to  geometrical  optics. 
It  contains  the  main  correction  term,  which  is  proportional  to  k   ,  corres- 
ponding to  the  second-order  WKB  approximation.   This  term  has  also  been 
computed  by  Keller,  Lewis  and  Seckler  (l7),  directly  by  the  WKB  method. 
Eq.  (6.23)  agrees  with  their  result  . 

The  complete  expression  for  the  wave  function  in  this  region,  according 
to  the  above  results,  is 
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^-W6  1/3/'   ^2  V/'^     exp[ik(r2-  a^)^/^: 

1         r        (2n)^/^       ^  Vr  -   aV  (kr   sin  9)^/^ 


X 


1 

"2 


[Ai'(-x    )1 
n 


I  ^    [Ai'(-x    )]'^         ^  ^ 

I  n  ^      n  m=l  n 

X      [exp(iX.^T^  +  i^)    +  exp(iX^5^-  i|)]L  (6.2l|) 

where  ih     and  ?//  are  given  by  (6.I3)  and  (6.23),  respectively.   From  the  point 
of  view  of  a  strict  asymptotic  expansion,  the  residue  series  should  not  appear 
in  this  expression,  since  it  is  exponentially  small  as  compared  with  the 
other  two  terms,  and  even  with  respect  to  higher-order  correction  terms  not 
taken  into  account  in  (6.23).   However,  from  a  physical  point  of  view,  this 
term  is  meaningful,  since  it  represents  the  continuation  of  the  surface 
waves  (5.6)  that  were  found  in  the  shadow  region.   Note  also  that  the  residue 
series  in  S  becomes  significant  near  the  shadow  boundary,  when  condition 
{6.k)   is  no  longer  satisfied. 


B.   The  Backward  Half-Space 

The  above  treatment  is  no  longer  valid  for  0  >  n/2.   Accor-^'ng  to  (^.8) 
and  ih.9),   the  integrand  of  (6.5)  blows  up  exponentially  for  |x|  ->  <»  in  the 
region  between  the  imaginary  axis  and  the  curve  t]  ->  -Tt/2  in  Fig.  6.   Thus, 
the  transformation  that  led  to  (6.6)  can  no  longer  be  performed.   Furthermore, 

the  integrand  of  (6.6)  no  longer  has  a  saddle  point  on  the  real  axis  between 

(2) 
P  and  p;  rather,  it  is  the  other  term  in  (6.5),  containing  }i\      (p),  that  has 

A 
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such  a  saddle  point. 

It  is  therefore  necessary  to  modify  the  above  procedure  in  the  whole 
backward  half -space.   For  this  purpose,  we  shall  again  start  from  (6.5),  but 
now  taking  the  path  of  integration  symmetrically  about  the  origin,  from  -oo  +  ie 
to  +00  -  ie  (note  that  there  are  no  poles  on  the  positive  real  axis),  and  we 
make  the  substitution  A  ^  -X  in  the  integral  from  0  to  +0°  -  ie.   Taking 
into  account  (2.15)  and  the  relation 

Q^^^cos  e)  -  e'^^''^^^^  (cos  0)   =  -e'^^'^an  nX  P   (-cos  9)  ,    (6.25) 

which  follows   from   (2.1ij-)  and   (C.5),   we   find 

-00  +  ie 
tp^-^'(T,Q)      =     2      /     f(X,3,p)P        (-cos    e)e'^^'^/^X  tan   rtX  dX    .  (6.26) 

o  ^-* 

The  behaviour  of  the  integrand  in  region  3  of  Fig.  h   differs  from  (4.7) 
(with  m  =  0)  only  by  the  replacement 

exp[iX(|  -e)]  -.       exp[-iX(|^-e)J  , 

so  that,    for   X  ^    aoo  (cf.    (2.26)),    it  behaves   like 

exp[-|x|(Ti   +  0   -   n)]    , 

where  r\   is  given  by  (2.27).   It  follows  that  the  path  of  integration  in  (6.26) 
may  be  deformed  into  the  path  from  0  to  ooo  shown  in  Fig.  3,   with 

n  >  JT  -  ©  .  (6.27) 
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Thus, 


^^^^(r,e; 


,  jt  ^2  -±n/h 


a  00 


Hf'(p) 


X  p   (-COS  e)e  ^^'^'^tan(rtX)XdX 


(6.28) 


In  particiilar,  for  0  >  7i/2,   we  may  subject  r\   to  the  additional 
restriction  (2.31).   Under  these  conditions,  (2.35)  applies  and  (6.28) 
becomes 


,{^\      rs\  IP  cos  9  ,  ,, 

t//^   (r,0)   =  e         +  f     , 
^  o  r 


(6.29) 


where 


^. 


/^   r   -77^ h(^\p)P^  i(-cos  0) 

o     X   ^ 


X  e'^^''/^  tan(rtX)X  dX 


(6.30) 


We  may  now  apply  again  (6.25)  and  the  converse  of  the  transformation 
that  led  from  (6.5)  to  (6.26).   The  result  is 


'^r      ^2p 


V  1  .)  /  o  \    X        X-  2 


^     H^^'O) 
o  «>     X 


(0  >  rt/2)  , 


(6.31) 


where  the  path  of  integration  is  that  of  Fig.  3  taken  in  the  opposite  sense. 


This  integral  may  be  evaluated  by  the  saddle-point  method.   There 
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is  now  a  single  saddle  point  on  the  real  axis,  at  A,  =  kp,  where  p  is 
again  given  by  (6.19).   In  fact,  the  integrand  is  identical  to  that  of 
(6.6),  the  only  difference  being  that  the  path  of  integration  now  goes 
over  only  one  saddle  point.   Thus,  the  result  of  the  saddle-point 
calculation  is  identical  to  (6.23). 

In  spite  of  the  fact  that  the  solution  in  the  backward  half-space 
is  just  the  continuation  of  the  solution  in  the  forward  half- space,  it 
does  not  seem  possible  to  extend  the  representation  {G.G)   to  0  >  n/2, 
or  to  extend  (6. 31)  to  0  <  jr/2.   The  reason  for  this  is  that  P   (x) 

becomes  singular  at  x  =  -1,  so  that  one  cannot  find  a  single  representation 

that  reinains  valid  both  for  0  =  0  and  for  0  =  rt.   One  needs  a  representation 

in  terms  of  P   (cos  0)  (cf.  (6.2))  near  0=0,  and  one  in  terms  of 
A.-2 

P  i(-cos  0)  (cf.  (6.28))  near  0  =  7t.   The  appearance  of  two  saddle 

X-2 

points  in  the  forward  half-space  is  also  related  to  the  diffraction 
effects  that  arise,  as  will  be  seen  in  Section  VII,  when  these  points 
approach  each  other. 

We  must  still  see  what  happens  to  the  residue  series  near  ©  =  n, 

f  2  )  f  2  ) 

because  both  the  term  ii^    '   containing  Q^  ^cos  0)  in  (6.2)  and  the  terms 

'//  with  m  >  1  in  (2.17)  become  singular  at  this  point.   However,  we  can 

write 
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00  +ie 


m=l 


f(\,P,p)Q^^''(cos   B)e^^''^^M\ 


-00  +ie 

00  CO   +i£ 

m=l  -00+  ie 


^4 


^^J(cos  e)+  Q^^|(cos  0) 

-    A,"  2^  A,"  2" 


X  exp[iTt\(2m  +  |)]XdX     =     2  \ 


oo  +ie 


nHI 


2  2^(-ir  J    f(x,i3,p: 

m=0  -00  +ie 


X 


^_i(cos   0)-   e  Q^_i(cos   0)      exp[irtX,(2m  +  ~)]ldX  ,         /g   ^2) 


and  employ  the  identity   (cf.    (C.5)  and   (C.6)) 


Q^^J(cos   0)-  e^^'^^Q^^^^cos   0)   =     -ie^"^     ^(-cos0) 

X.-'2  X~2^  X.~2^ 


(6.33) 


to  get 

00  00  oo 

'^^o^   ^    )K     =     aiV    (-if    rf(\,P,p)P     J-cos   0)exp[irt\(2ra-  \)]m\  . 
™=^  "^=^  -="  ^"'  (6.3I1) 

The  contour   of  integration  may  be  closed   in  the  upper  half-plane, 
leading  to  the  residue   series 


1  00 


^res^^^e)      =     2ne-^^/\§)'^    (-if^    Vn^xp[i«V'2m  + 1)] 


m=0  n=l 


.(1), 


X      H^     (p)P       ^(-cos   0; 
n  X  -2 


=      «e  (— ) 


2p'       ^    Vn 
n=l 


exp(ijrX  /2)        /-|_v 

cos   nX  4     (P)P,      J-^°^   ®)    ^ 


n  n 


X  -2 
n 


(6.35) 


/. 

T, 

^^/^ 

P(r,^)^ 

\ 

--M^^^T/ 

h\ 

~f^ 

\ 

»  7 

\ 

\    cos"'  7-y 
\\          /a 

v) 

/ 

^^ 

_^y 

T. 
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which  differs  from  (5.3)  only  by  the  substitution 


P   ^(cos  0)  ^  -i  exp(ijt\^)P   ^(-cos  9)  .  (6.36) 

A,  -"2  ■^  ""2 

n  n 


If 


rt  -  0  »  P'^  ,  (6.37) 


we  may  employ  the  expansion  (C.8)  for  P    (-cos  0)  and  we  find  that 

n 
(6.35)  is  equivalent  to  the  residue  series  appearing  in  (6.2k),    so  that 

{6.2k)   may  be  continued  to  0  >  n/2. 

Note  that 

T   =  2n  +  -  -   cos"  (a/r)  -  0 
1  2        \    /     ' 

no  longer  corresponds  to  one  complete  encirclement  of  the  sphere,  but 
rather  to  the  smallest  angle  described  by  a  sirrface  wave  excited  at 
T-,  before  leaving  the  surface  (Fig.  8). 

If,  instead  of  (6.37),  we  have  «  -  0  ~  P   ,  we  must  employ 
(C.9)j  which  leads  to 

m=0       n    L  i^  n'-J 
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This  expression  snould  be  compared  with  (5.II).   The  same  focussing 
effect  already  discussed  there  leads  to  an  enhancement  of  the  radiation 
from  the  surface  waves  in  the  backward  direction. 

The  approximations  employed  in  the  present  Section  fail  when 
the  correction  terms  appearing  in  (6.23)  become  large.   This  happens 
near  the  forward  and  the  backward  directions.   The  behaviour  near  the 
backward  direction  will  be  discussed  in  Section  IX,  in  connection  with 
the  scattering  amplitude.   The  bejiaviour  near  the  forward  direction  and 
the  corresponding  diffraction  effects  will  be  investigated  in  the  next 
two  sections. 
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VII.   Diffraction  effects  in  the  near  region 

Ac   The  Normal  Derivative  on  the  Siarface 

According  to  Huygens '  Principle,  the  wave  function  at  any  point  in 
space  can  be  expressed  in  terms  of  its  normal  derivative  on  the  surface  of 
the  sphere  (where  the  wave  function  vanishes). 

The  analogue  of  Kirchhoff's  approximation  in  classical  diffraction 
theory  would  he  to  replace  the  exact  values  of  the  normal  derivative  by  the 
geometrical  optics  approximation: 

X  =  0         (0  <  rt/2)  ,  (7.1) 


^  -  I 


B  /  ikr  cos  0, 


o-        ^   iP  cos  0  /^  .   JT, 

=  2i  cos  0  e         [9  >  -, 
r=a 


(7.2) 


where  we  have  introduced  the  notation 

The  expression  (7.1)  corresponds  to  the  geometrical  shadow  region  and  (7.2) 
to  the  geometrically  lit  region  on  the  surface.   The  factor  2  in  (7.2)  arises 
from  joining  the  contributions  of  the  incident  and  geometrically  reflected 
waves. 

Substituting  (7.1)  and  (7.2)  in  the  Huygens -Kirchhoff  integral,  one 
can  easily  derive  the  corresponding  approximation  for  the  scattering  amplitude 
(l8).   It  is  fo\md  to  be  the  sim  of  two  terms,  one  of  which  corresponds  to 
the  geometrically  reflected  wave,  while  the  other  one  corresponds  to  the 
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diffracted  wave  in  the  classical  theory  of  diffraction  by  a  circular  disc  of 
radius  a.   The  latter  term,  which  is  also  known  as  the  "shadow-forming  wave", 
depends  only  on  the  shadow  contour  of  the  obstacle,  so  that  it  is  the  same  for  a 
sphere  or  for  a  disc. 

Since  we  are  later  going  to  compare  ovor   resiilts  with  classical 
diffraction  theory,  it  is  of  interest  to  discuss  the  accuracy  of  Kirchhoff's 
approximation  by  evaluating  X.   The  only  significant  contribution  arises 
from  (6.5): 


X 


(1) 


00  +ie  Q^  ^(cos  e) 
-00  +ie     X     ^ 


il-^) 


where  we  have  employed  (2.8)  and  the  Wronskian  relation 


w[h[^)(P),  h[2)(3)]  .  -(4iA3)  . 


(7.5) 


For  e  >  rt/2,  the  main  contribution  to  (7.^)  arises  from  a  saddle  point 
on  the  real  axis,  to  the  left  of  A,  =  P,  where,  according  to  (A.16)  and  (C.7), 


,,eiW2^^)  ^^^^  9)/h(^V)  z     i(Vsln  9)^(3^-  X^  ^ 


x4 


exp 


-i(3^-  x.^)2  +  ix(|  -e  +cos"^(Vp) 


(7.6) 


so  that  the  saddle-point  is  at 


cos"^(>:/p) 


0  -  J  -*  X  =  P  sin 


(7.7: 


The  saddle  point  evaluation  gives 


1+2 


X^^^  Z     2i  cos  e  e^P  ^°^  ®     G  -  f  »  3"3^  ^   ^^_Q^ 


in  agreement  with  the  geometrical  optics  approximation  (7-2). 

The  condition  in  parentheses  arises  from  the  fact  that  (A.I6)  is 
valid  only  for  P  -  X  »  P  '  .   As  0  approaches  n/2,   the  saddle  point  (7.7) 
moves  towards  A,  =  P  and^  if  |0  -  -|   <  P    ,  we  have  to  employ  the 
approximation  (A. 17 ) : 

h[^\p)  ;=  2(2/p)3e-^"/3Ai(xe2^'^/3)       (|x-3|  <  P^  J ,      (7-9) 

where 

1/3 
X  =  (2/p)   (x-p).  (T-10) 

Since  the  main  contribution  to  the  integral   in  this   case  arises   from 

|x|   5  1 J  we  may  extend  the  range  of  integration  where   (7-9)   is   employed  to 

inf inity,   with  the   following  resixLt : 

,^,  -ijT/6  exp[ip(f  -  9)1  /  -l/3\ 

^°'    =     -T^i   -         1/3         ^'-'      {\e-l\  <,'-),    (7.11) 


where 


and 


=    (P/2)^/3   /|.e^  ^^_^2) 


ITX  p  -ITX 

e  n         ,       /  e 


':^(t)      =  -    '      ■    /  -  dx     +  — ^— r-77-     dx 

^  J        Ai(xe^"/3)  ^        Ai(xe-^^ 


o 


00  exp(2ijr/3)  -in/6 

ijr/3       r  exp(e        '^^    rw) 

j  Ai(w)  ^^  •  (7.13) 

«>  exp(-i7t/3) 
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Finally,  if  -  -  0  »  P    ,   we  may  evaluate  (T-^)  "by  closing  the 
path  of  integration  in  the  upper  half -plane,  where  the  integrand  is 
exponentially  decreasing  at  infinity.   This  leads  to  a  residue  series  at 
the  poles  of  [h^   (3)]   •   The  terms  in  the  series  are  then  rapidly 
decreasing  and  the  main  contriliution  arises  from  the  poles  (3.5): 


(sm  0)'^   '—'       '   n 


This  corresponds  to  the  surface  waves- associated  with  T  in  (5-6).   The  normal 
derivative  is  exponentially  damped  in  the  shadow  region  p  "  ©  »  3    , 
in  good  agreement  with  (7.1). 

Thus,  the  only  domain  where  Kirchhoff's  approximation  (7.1)  -  (7-2) 

fails  to  be  acciu-ate  is  the  penumbra  region 

_1 
|0  -  |I  <  P  3  (7.15) 

I 
where  the  normal  derivative  is  given  by  (7.11).   The  fimction  '^^(t),  which 

gives  the  transition  from  light  to  shadow,  was  introduced  by  Fock  (3_). 

There  remains  to  be  shown  that  (7.11)  goes  over  smoothly  into  (7.8) 
or  (7.1^^-)  for  |©  -  5|   »  P~"^'  .   The  asymptotic  behavioiir  of  the  integrand 
of  (7.13)  in  the  w-plane  follows  from  Appendix  D.   It  is  fo\md  that  the  path 
of  integration,  represented  by  the  straight  line  D  in  Fig.  S,   may  be  deformed 
at  will,  provided  that  it  begins  and  ends  at  infinity  outside  of  the  shaded 
sector. 

If  T  »  1,  we  may  evaluate  (7.13)  by  moving  the  path  of  integration  to 
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infinity  in  the  left  half -plane,  which  leads  to  a  residue  series  at  the  poles  -x  ; 

/  -iTi/6        ■. 

■    la.  'S""'  exp(-e   '  tx  ) 

^(T)   =  2^.-^^l^2_,  Ai-(-x  )   "  (t»1).     (7.16) 

n 
On  the  other  hand,  if  t  «  -1,  we  may  evaluate  (7.13)  liy  the  saddle-point 

method.   For  this  purpose,  D  is  deformed  to  the  right  over  the  region  where 

|w|  »  1,  so  that  the  expansion  (D.  ij-)  can  be  employed.   The  saddle  point  is 

located  at 

and  the  corresponding  path  of  integration  C  is  shown  in  Fig.  9.   The  result  of 
the  saddle -point  evaluation  is 

:^{^)    %    i^jTe'^''/3|^j  g^p(.|^|3/3^  (t  «  -1)     (7.18) 

Substituting  (7.l6)  in  (7-11),  we  get  (7.1^).   Substituting  (7.18) 
in  (7.11),  we  find 


fl)        (|-®) 
X^-^^  s  2i  -^ T  exp.(ip 


O  /  .   .^2 


;sin  e)' 


(|-e)-i(f-e)3 


(7.19) 


which  agrees  with  (7.8)  for  0  -  -  »  p  '   (but  still  not  too  large,  so  that  P  cos  9 
is  well  approximated  by  the  expansion  within  ciur-ly  brackets  in  (7.19)  ). 

Thus,  Fock's  function  indeed  interpolates  smoothly  between  the  shadow  and 
lit  regions  on  the  surface  of  the  sphere.   However,  it  cannot  be  employed  in 
the  lit  region  too  far  beyond  ©  -  "  ~  P    . 


Im  w 


Rew 


Fig.    9 
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B.   The  Neighborhood  of  the  Shadow  Boundary- 
Let  us  now  consider  the  behaviotir  of  the  wave  function  in  the 
neighbourhood  of  the  geometrical  shadow  boundary,   0  si  0   (cf.  ('^.2)), 
at  not  too  large  a  distance  from  the  sphere, 

r  «  p^/\  .  (7.20) 

Under  these  conditions  it  follows  from  (5.8)  and  (5.9)  that 

|exp(iX^&^)|   «  1  ,  (7.21) 

so  that  we  need  only  be  concerned  with  t   ,  which  is  given  by  (6.6). 

We  have  seen  that  (6.6)   has  saddle  points  at  X  = P  sin  0  and 
\  =  kp,  where  p  is  given  by  (6.19).   The  neighbourhood  of  the  shadow 
boundary  is  characterized  by  the  fact  that  these  two  saddle  points 
approach  each  other,  moving  toward  the  point  X  =  P.   The  main  contribution 
to  the  integral  then  arises  from  the  neighbourhood  of  X  =  P,  so  that  it 
is  convenient  to  take  the  path  of  integration  through  this  point  and  to 
split  it  in  the  following  form: 


* 


(1)  ,     .   (^)i  e-W4  r  [     hJH   „(l'(p„(l'  (cos  e)  e'^"/2; 
,  ,      r   iii!i  „(l)(p)  «(l)(eos  e)eiW2^, 

J  h[^'(3)  ^  x-i 

-  J  hU)(p,  ,(i)j,o3  e,e-/^x.xj  -    C^*l','^*l;'. 

P  X-  2  -^ 

(7.22) 


hG 


where  a  »  and  a^eo   denote  the  beginning  and  end  points  of  V,   as  illustrated 
in  Fig,  G,    and  we  have  used  the  identity 

h[^^(P)/h[^V)  =  2j^(3)/h[^^(P)  -  1  .  (7.23) 

In  ti-5  J  we  may  employ  the  same  approximations  that  led  to  (6.9). 
Restricting  ourselves  to  the  main  term  of  (6.11),  we  find 

1"  .  ^T.24) 

i\,^       ss      (— : — -)^     e        '         /    (sin  w  cos  w)^expfip  [cos  w  +  (w  -0)sin  w] )  dw 

03  2n   sm  0  J 

e 


Above  the  shadow  boundary,  the  integrand  still  has  a  saddle  point 

at  w  =  9,  as  in  (6.10).   In  any  case,  for  G  ~  0  ,  the  main  contribution 

comes  from  the  neighbourhood  of  the  lower  limit  of  integration,  so  that 

we  may  expand  the  integrand  around  w  =  0  and  extend  the  corresponding 

range  of  integration  to  infinity,  with  the  following  result: 

i_  00 

0 
o 

where  we  have  substituted  p  cos  0  =  kz. 

The  condition  for  the  validity  of  this  approximation  is  that 
higher-order  terms  in  the  expansion  of  the  exponent  in  (7.2^)  shall  be 
negligible  in  the  relevant  portion  of  the  domain  of  integration.   This 
leads  to  the  following  conditions: 

kz  »  p2/3^  (7.26) 

|0  -  0^1  «  p"^/3^  (7.27) 


where 


and 


hi 


The  above  resiilt  may  be  rewritten  as  follows : 


^(1)   ^  ^ikz-i.A^^^^)  -F(-v)]A/^  (7.28; 


V      =      (kz/jt)2  (e  -  9  )  (7.29) 


F(v)  =   /  exp(ijtT^/2)dT  (7.30) 

o 

is  the  Fresnel  integral.   Since 

F(co)   =  e^'^/Vv^  ,  (7.31) 

(7.28)  becomes 

r-,  \  1       ■^  ikz-in/ij- 

^^^)  .     |e^^^  -^—r—     F(-v)  .  (7.32) 

•^  y2 

This  is  analogous  to  the  classical  Fresnel  diffraction  pattern  of 
a  straight  edge  (19).   On  the  geometrical  shadow  boimdary,  where  v  =  0 
(cf.  (7.29)),  it  woiild  give  one  half  the  amplitude  of  the  incident  wave, 
which  is  a  well-known  result. 

Now  let  us  consider  the  two  remaining  terms  of  (7.22).   Since  the 
main  contribution  to  the  integrals  comes  from  the  neighbourhood  of  X  =  P, 
we  may  employ  the  approximations  (A. 17)  and  (A.I8),  which  give 

h[2)0)/h[1)(P)  ==   e2W3Ai(xe-2i"/3)/Ai(xe2i«/3)  (X<P),   (7.33) 


2JJP)/h[^^(P)   ~  e^"/\i(x)/Ai(xe2^''/3)  (X>3),   (7.3^^) 


l(l)O)   ^  e^"/3Ai(xWA.^..2in/3. 

A 

where  x  is  defined  by  (7.10). 
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For  the  remaining  factors  of  the  integrand,  we  may  employ  the 
same  approximations  that  led  to  {1.2^-),    so  that  we  get,  with  the 
substitution  (6.8), 

0 


^01^    ^^02^ 


■c 


2n  sin  9 


2      -±n/k 
)      e        ' 


2i:t/3       r'    .  i  Ai(xe-^^"/3) 

'         '    (sm  w  cos  w)     — '   •\,. — r- —   . 

2iit/3. 


X  exp  <  ip[ cos  w  + (w  - 0)sin  wj    >  dw 


(sm 


Ai(xe 


.2  Ai(x) 

W    cos    w) o  •      /■:."  ' 

„  .  ,      2ijt/3  \ 
Ai(xe        '     ) 


X 


exp  i/ip[cos  w  +  (w-  0)sin  w]    >  dw 


:7.35; 


In  the  neighboiirhood  of  0  =  0  ,  we  may  again  employ  the  expansion 

around  w  =  0  that  led  to  (7.25),  with  the  following  result: 

0 
o 


01  02 


2)T  " 


■  (ol )     e        '       e         <  e 


exp[ikz(w-  0)   /2] 


Aifxe-^^"/3)      ^        ^ 
X       — '   ■  ,-,. — 7- —     dw     +     e 

Ai(xe        '     ) 


in/3      r         r-i     ^        Q^2,         Ai(x) 

'         /   exp[ikz(w-0)   /2J '— ^ 


Ai   xe 


2iit73^ 


dw 


y,  (7.36) 


where 


x      = 


2  1/3 
(o")        p(sin  w  -  sin   0 


2  1/3  /'^"^^^  ,''^'  ®o^ 

2(r)       p  cos  i       g  °)   sin  (-- r— )  .      (7-37) 


The  discussion  of  the  asj^.ptotic  "behavioiir   of  the  Airy  function  given 

in  Appendix  D  shows  that  values   of    |x|    »  1  do  not   contribute  much  to  the 

1  /^  2./^ 

integrals   in   (7-36).      On  the   other  hand,    p/|3~'-   >  p   '-^  »  1,    so  that. 
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according  to  (7.37):,  only  small  values  of  |w  -  0  |  give  a  significant 
contribution.   Under  these  conditions,  we  may  replace  (7-37)  by 

X     ^  (2/(3)^/3p  ^o3  0  /^_Q J  ^  (2/p)^\z(w-e^)  ,       (7.38: 


where  the  last  approximation  follows  from  0  s;  0. 

Making  the  change  of  variable  (7-38)  in  (7-36),  we  find 


o 

,(1)        ,(1)  -i     -irt/U       ikz        ,.      2-    , 

^01  ^02       ~     ■'^       ^        '      ~ne       exp(iJtv   /2 ) 


2i7r/3    r        ft     ^-22, 
/    exp(i§x  +  iT^  X    , 


(7.39) 


Ai(xe"-^"/^)      ^      ^     i«/3    r  f.,      ^.22.  Aifx) 

^     o.-^/?  dx  +  e     '-"    /      exp  i|x  +  n^   x    ) ^TtTT  '^^ 

Ai(xe2^''/3)  ^  Ai(xe2^"/3) 


where   v   is  given  by   ('J. 29),   and 

i      =      (p/2)^/3(Q   .   Qj^  ^      ^      (p/2)^/V(2kz)^    .  (l.kO) 

According  to  (7.26)  and  (7-27),  we  have  \  i\    «  1,  f]   «  1,  so  that,  in 
the  significant  range  of  the  domain  of  integration,  the  exponential  function  in 
the  integrands  of  (7.39)  may  be  replaced  by  unity,  with  the  following  result: 

^01^  ^^02^  ~  -Ce"^^/^  P^/3(2nkz)"2exp(ikz  +  im.^/2)  ,  (7.i^l) 

where 

,1/3,  ,  ,21,/3  J    Ai(xe-^;;f '  ,,.,W3  r       -^iljV...   (7.^) 

The  constant  C  has  been  evaluated  by  Rubinow  and  Wu  (^)  who  found 
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0.99615  e^"/3 


(7.U3: 


It  follows  from  (7.28)  and  (7.^)  that 


t 


(1) 


e       ikz 

e 


v^ 


F(oo)-F(-v)-2l/3   exELi2v_Z2_ 


(7.Ulf) 


The  conditions  for  the  validity  of  this  result,  besides  (7.27),  are  (7-20) 
and  (7.26),  which  may  be  combined  into  the  following  condition: 


p  ^'^  a   «  z  «  p-^' 3  a 


(7.^5) 


It  is  readily  seen  that,  for  |v|  »  1,  either  in  the  lit  or  in  the 
shadow  region,  the  last  term  of  (7.^^)  is  a  small  correction,  of  the  order 
of  [|[,  to  the  amplitude  of  the  diffracted  wave. 

Thus,  if  |0  -  0  I  «  p   '-',  in  the  domain  (7-^5),  the  transition  from 
light  to  shadow  is  described  by  an  angular  Fresnel  diffraction  pattern  very 
similar  to  the  classical  one  for  a  straight  edge.   The  effects  of  the  curvature 
of  the  sphere  come  in  through  small  correction  terms,  of  the  order  of  \i\. 
If  1 0  -  0  I  »  p    ,  we  go  over  either  to  the  lit  region,  where  (6.23)  is 
valid,  or  to  the  shadow  region,  where  (5-6)  is  valid  (cf.  (5. 10)). 

Finally,  let  us  consider  the  immediate  neighbourhood  of  the  shadow 
boundary,  where  (cf.  (7.29)) 


«  1  . 


(7.it6) 


We  may  then  approximate  the  Fresnel  integral  by  the  first  term  in  its  power 
series  expansion. 
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p  cos  e(e-e  )     D(sin  e  -  sin  0  ) 

-F(-v)  «  V   =  T ~ T ,     (7.^7) 

(rtkz)2  (jtkz)2 

and  the  exponential  in  the  last  term  of  (7-^^)  may  be  replaced  by  unity. 
Taking  into  accoimt  (7.31  )j  we  get 


,(1)        1  ikz 
*o       ~F 


/^  2  \  ^      -±n/k  ,  ,        .      ^  no-2/3    N 

1  +   [~r~)       e        '      k(r  sm  9-  a-  Cp     'a) 


(7.^; 


Defining  the  shadow  boundary  by  the  condition  that  \ii\    =  —  on  it,  we  find 
from  (7-^)  that  the  shadow  boimdary  no  longer  lies  at  r  sin  0  =  a,  as  in  (7.32), 
but  rather  at 

r  sin  e  =  a  +  s  ,  (7.^9) 

where 

s  =  a(Re  C  +  Im  O/p^^^  =  1. 36077a/P^/^  .  (7-50) 

This  is  identical  to  the  result  foimd  by  Rubinow  and  Keller  (21^)  for  a 
circular  cylinder  and  by  S.  0.  Rice  (22)  for  a  parabolic  cylinder,  thus 
confirming,  in  the  present  example,  Rubinow  and  Keller's  conjecture  that  the 
result  is  true  also  for  three-dimensional  obstacles. 

There  remains  to  examine  the  consistency  of  the  various  approximations 
leading  to  (7.^).  According  to  (7.^7),  (7.^9)  and  (7-50),  the  order  of  magnitude 
of  V  at  the  (shifted)  shadow  boundary  is  given  by 

ill 

(nkz)2      (jtkz)' 


_±S^^     ~  -§-— X  (7.51) 

\2  f^-U'7\2 


According  to  (7.26),  this  satisfies  condition  (7.^)- 
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It  is  readily  verified  that,  in  the  domain  defined  by  (7.27)  and 
(7.ij-5),  the  neglected  terms  in  (6.11)  and  in  the  approximations  that  were 
made  in  connection  with  (7-35)  give  contributions  of  a  higher  order  of 
magnitude. 

VIII.   The  Fresnel-Lommel  region 

A.   Basic  Approximations 

According  to  classical  diffraction  theory,  the  Fresnel  region  is  the 
domain,  in  the  neighbourhood  of  the  shadow  boundary,  viewed  from  which  the 
obstacle  contains  a  large  number  of  Fresnel  zones,  i.e.  r  «  a  /X,  where  X 
is  the  wavelength.   We  shall  now  consider  the  behaviour  of  the  wave  function 
in  this  region,  at  distances  larger  than  those  allowed  by  (7.20),  i.e. 

P^/\  <  r  «  Pa  .  (8.1) 

At  the  same  time,  we  shall  stay  within  the  geometrical  shadow,  or  not  too  far 
outside: 

e  <  9   -  a/r.  (8.2) 

~   o 

We  shall  call  this  the  Fresnel-Lommel  region,  because,  as  will  be  seen 
later,  the  wave  function  throughout  most  of  this  region  can  be  approximated  by 
Lommel's  classical  solution  to  the  problem  of  diffraction  by  a  circular  disc (23). 

According  to  (5.8),  (5.9),  (8.1)  and  (8.2),  condition  (7-21)  is  no  longer 
satisfied  in  this  region,  but  rather  |exp(iX  &  )|  >  1.   Thus,  the  integral 
containing  Q,  i(cos  9)  in  (6.2)  can  no  longer  be  expressed  as  a  residue  series 

X-2 

and  we  must  consider,  in  the  place  of  (7.22), 


1 
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•a^oo   h;^  ^'p) 


J^—   H[^^(p)P^_i(cos  e)e^^/2^^ 


"'  /  MTZ:     4'^^^)  Vi^-^  e)eiW2^,         ^g^3) 
3   \  ^^^ 


O^oo 


3 


HJ^^V)  Pj^.iCcos  Q)e^^/^MX 


=     toi  +  to2  +  *03 


The  third  term  of  (8.3),  like  that  of  (7.22),  depends  on  the  scatterer 
only  through  its  radius  appearing  in  the  lower  limit  of  integration.   Roughly, 
it  represents  the  effect  of  cutting  off  from  the  incident  wave  all  the  rays 
that  meet  the  sphere,  so  that  the  sphere  hehaves,  in  this  respect,  like  an 
opaque  disc  of  radius  a.   We  shall  see  that  this  term  gives  rise  to  the 
classical  diffraction  pattern  of  a  circular  disc.   It  corresponds  to  the 
"shadow-forming  wave"  mentioned  in  Section  VII. 

The  main  contribution  to  the  integrals  in  (8.3)  still  arises  from  the 
neighbourhood  of  X.  =  P,  so  that  H;   (p)  may  be  replaced  by  the  expansion  (A.l6). 

A. 

On  the  other  hand,  since  we  want  to  consider  both  the  behaviour  for  P©  <  1 

and  for  p0  »  1,  we  shall  employ  the  uniform  asymptotic  expansion  (C.ll'l  of 

P  x(cos  0).   Finally  in  y\f^^    and  i'^,   we  may  again  employ  the  approximations 
A,-  2  U±       Od 

(7.33)  and  (7-3^).   The  results  are 


5U  - 


^01  ^  ^02 


D  \sin   9/ 


/3      r      Ailxefi!/ll 


2        2, 


a  ,00 


Ai(xe 


21^75,      ^"^^^^^    "  ^  ^ 


CTgOo 


+  iX  sin"      (x/p)]j   (A,e)Xd\  +  e^ 


rt/3      r  Aifx)  r  ,   2       2.i 


+     i  \  s  in      (  A,/p  ^ 


J^(A,e)X.dX.>     , 


(8A) 


where  x  is  given  by   (7.10)   and 


0^00 


'03 


-ip(e/sin  9)^      /   exp^ip[(l-   r    f  +  t   sin      t]  V  J   (peT)(l-T    )      rdr  , 

a/r  I  J      °  (8.5) 


1 


where  we  have  made  X  =  pr. 


Since  the  main  contribution  to  the  integral  in  (8.5)  comes  from  the 
neighbo-urhood  of  the  lower  limit  t  ~  a/r  ~  P  '  «  1,  we  may  expand  the 


integrand  in  powers  of  r,   keeping  only  the  main  terms: 


'03 


-ip(e/sin  e)^  e 


2  „iP 


exp(ipT  /2)J  (p0T)TdT  , 


(8.6) 


a/r 


where  the  upper  limit  has  been  replaced  by  00. 


According  to  Appendix  E,  this  integral  may  be  expressed  in  terms  of 
Lommel  functions.   In  fact,  it  follows  from  (E.5)  that 


]^       s  (o/sin  0)^exp(ip  +  ik  — 


03 


2r 


y^i~-,^Q)   +  iv^(— ,pe) 


(8.7) 


where  V  and  V.,  are  Lommel  functions  of  orders  zero  and  one. 
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Similarly,    in   (8.4),   we  may  approximate 


(p      -    \   )^   +   X  sin'    (X./p) 


+    (A.72p)    . 


(8.8; 


Taking  x  as  new  variable  of  integration,  we  find 

X-/2P   =  ka2/2r  + '0/2  )^/\x/r  +  iii/2f^\^/2p    .  (8.9) 

The  main  contribution  to  the  integrals  comes  from  ix|  ^1,  because  of  the  Airy 
functions.   Thus,  according  to  (8.1),  the  last  term  of  (8.9)  is  negligible,  and 
we  find 


.  2i«/3  r  Aifxe-^""/3^ 


+  e  ' 


Ai(xe 


Ai(x) 
Ai(xe   '  ) 


i(|)^/3j. 


B  1/3  -1 
3G  +  (|)    0x 


dx 


(8.10) 


i(|)^^3  J  X 


136  +  (|)^/3g^ 


dx, 


Putting  together  (8.7)  and  (8.10),  we  finally  get 


9 


.sin  0 


exp  (ikr  +  ik  ~  j   f(s,t,u,v)  , 


(8.ii: 


where 


f(s,t,u,v)   =  L(u,v)  +  isF(s,t,v)  , 


L(u,v)   =  V^(u,v)  +  iV^(u,v)  , 


(8.12) 
(8.13) 


o 
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-2in/3. 


.2i«/3      r   Aifxe;-^"^^)      ^isx 


F(s,t,v)      =     e-^^^^^      /     -   ^      2i^/3   ^     e---  Jjv   +  tx)dx 


(8.U) 


and 


J     Ai(x7^3) 


s   =   0/2)^/\/r,  t   =   0/2)1/3  0^ 


u    =    ka   /r      =     Pa/r  ^      v  =  P0  =  tu/s 


(8.15) 


Conditions  (8.1)  and  (8.2)  are  equivalent  to 

s  <  1,    t  <  1,   u  »  1  .  (8.16) 

According  to  classical  diffraction  theory,  the  wave  function  in  the 
P^esnel  region  due  to  the  diffraction  of  a  plane  wave  by  a  circular  disc  of 
radius  a  is  (23 ) 


^cl 


exp  (ikr  +  ik  ^j      V^  (^,    P  sin  9J  +  ^\(^>   P  ^i"  9\   ]  ,   (8.1?; 


in  the  approximation  where  sin  0  ~  0.   In  this  approximation,  (8.17)  coincides 
with  (8.7). 

Thus,  in  (8.12),  L(u,v)  represents  Lommel  s  approximation,  while  F(s,t,v) 
is  a  correction  to  classical  diffraction  theory  of  the  same  type  as  the  Fock 
terms  discussed  in  Section  VII. 
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B.   Behavioirr  on  the  Axis  (Poisson  Spot) 

For  0  =  0  (r  =  z),  we  have,  according  to  (E.3)^ 


so  that 


where 


L(u,0)   =  1, 


\|f  ~  exp  f  ikz  +  ik  ~  )  f(s)  , 


f(s  )   =  1  +  is 


■^2i./3  f     Aifxe^A   ^-isx^^ 
^o   Ai(xe-«/3) 


+  e   '  -^ 


Aifx) 

„  .  ,   2i3t/3  , 
Ai(xe   '  . 


isx^ 
e   dx 


(8.18) 
(8.19) 


(8.20) 


In  Lommel's  approximation,  the  second  term  would  be  absent,  so 
that  we  would  have  |t  |  =1?  i.e.  the  intensity  along  the  axis  would  be 
identical  to  that  of  the  incident  wave.   This  corresponds  to  the 
well-known  Poisson  spot. 

According  to  Appendix  D,  Ai(^e  '     )/kl{^e       '     )   goes  to  zero  like 


exp 


\    U|3/2sin(3-|72 


when  C  =  I  ^|exp(i'45)  goes  to  infinity  in  the  sector  kn/Z  <  'J?  <  2-n   so  that 
the  path  of  integration  in  the  first  integral  of  (8.20)  may  be  rotated 
by  jt/3  into  the  fourth  quadrant,  with  the  following  result: 

00 

_iji/3" 


f(s)   =  1  +  ise 


Ai(x)  ,  _ 

Ai (xe' 


ii/JT  -p(-i<=-"^'-to 


Ai(x) f.       .^' 

A.  /o   exp(isx)dx 

AKxe^^'^) 


(8.21) 
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/  ^  3/2\ 
Both  integrands  now  tend  to  zero  like  exp  {  -  rx  '  1   for  x  -»  ». 

At  distances  z  «  3   a  (s  »  1),  the  residue  series  representation 
(5.11)  should  converge  well,  so  that  the  wave  function  must  'become 
exponentially  small.   To  recover  this  result  from  (8.21)  let  us  note 
that,  by  partial  integration,  (8.21)  becomes 


f(3)   =  e___   /  exp(-ie-  /  sx)     ^  e^L  T  __exE(i|x)     ^^   ^g^^^) 
^         2«    j    [Ai(xe-2^"/3)]2       2«   J   [M(xe2^"/3)]2 
o  o 

where  we  have  employed  the  Wronskian  relation  (D.2).   This  result  may  be 

rewritten  as 

f(s)   =  ^^  f      — ""^^-^H   2^^-         (^-23) 

where  P  is  the  path  shown  in  Fig.  10,  going  from  ooe   '   to  00. 

If  s  »  1,  the  integral  can  be  evaluated  by  reducing  it  to  a  series 

of  residues  at  the  poles  ^  =  e""""'  x  ,  where  -x  is  the  nth  zero  of  the 

■^      n        n         n 

Airy  function.   The  result  is 

/ .    ijt/3        \ 
.     /r       \. — 1     exp(ie  sx    ) 

f(s)      =     e-^'^/^s      > f-'  (8.2U) 

^-'        [Ai'(-x    )Y 
n  L        ^      n  -^ 

Substituting  this   in   (8.19),   we  find 

.    /.  /  2  \  ,,  .  1/3  ^     V^  exp[ie^"/3      (p/2)l/3a/z] 

(z  «  3^/3a).  (8.25) 

It  may  readily  be  verified  that  this  result   coincides  with  the  term 


o 


Fig.    10 
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m  =  0  of  (5. 11) J  with  z  »  a.   Thus,  the  wave  ftinction  along  the  axis 

1/3 
is  exponentially  small  for  z  «  P  '  a. 

1/3 
On  the  other  hand,  for  z  »  P  '  a  (s  «  1),  we  can  expand  f(s) 

in  a  power  series, 

f(s)   =  1  +  is  \  i"M^s"/n:  (8.26) 

n=G 

where,  according  to  (8.21), 

00  00 

M   =     e^^/3   T-^iM^  x"dx  +  ei^2n+l)./3  f   ^Ail^.  ^n^^ 

J      Ai(xe2^"/3)   •  J      Ai(xe-2W3) 

o  o    ^ 

(8.27) 

The  first  few  coefficients  M  have  "been  computed  by  Wu  (24).   In 
particular, 

M^   =  2^/3^   ^  1.2551e^''/3,      M^   =  0.5323e^^''/3,      (8.28) 

where   C  is   given  by   (7.^3). 

Substituting    (8.26)   in   (8.19),   we  get 

/  2n    r  1/3  2/32  -, 

t^     ^  exp   f ikz  +  ik  I7  )    l  +  ±\{f)         I  -  \(P         \  +   --- 

(z  »  p^/\).  (8.29) 

1/3 
Thus,  for  z  »  p  '  a,  the  intensity  approaches  that  of  the  incident  wave. 

On  comparing  (8.29)  with  (8.25),  we  see  that  a  Poisson  spot  of 
intensity  comparable  to  that  of  the  incident  wave  begins  to  develop  at 
a  distance 

z  ~  pl/\  (8.30) 
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from  the  sphere.   This  is  in  contrast  to  the  case  of  a  circular  disc, 
where  the  Poisson  spot  begins  to  appear  at  a  distance  z  ~  a  from  the 
plane  of  the  disc  (5,  p.  103). 


C.   Behavioior  Away  from  the  Axis 

Let  us  consider  first  the  term  L(u,v),  which  corresponds  to 
the  classical  diffraction  pattern  (cf.  (8.13)).   The  hehaviour  of  this 
term  depends  on  the  parameter 

u/v  =  a/re  s;  e  /e  .  (8.31) 

In  the  lit  region,  9  »  9  ,  we  can  employ  the  expansion  (E.l), 


which  leads  to 
L(u,v)   =  exp 


1  /      V 


2  T  2 

u 


+  i  ^  J^(v)  +  ip    J^{v)    +  ... 

(u  «  v).    (8.32) 


The  corresponding  classical  wave  fixnction,  according  to  (8.11),  would 
be 


i' 


Q 


c.l  ~  Vsin  e 


exp 


ip(l-%) 


1  —  exp  (ip  +  ik  ^-)J^{^e)  +  ■■ 


exp(ip  cos  9)     (0»a/r).    (8.33^ 


so  that  if  .     approaches  the  incident  wave. 

Near  the  shadow  boundary,  at  9  =  a/r,  we  can  employ  (E.2),  which 
gives 


L(u,u; 
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so  that 


t 


e 


c^    2  Vsin  e 


exp  lip 


2\ 
•1  s 


)h-(-v>o(T^); 


(e  =  a/r). 


(8.35: 


The  factor  l/2  corresponds  to  the  classical  behavioijr  at  the  shadow  boundary 

1/3 
discussed  in  Section  VII.   For  r  »  P   a,  the  shift  of  the  shadow  boundary 

is  no  longer  given  by  (7-50):   it  increases  with  r  and  then  oscillates.   At 

distances  r  ~  3a,  which  mark  the  transition  to  the  Fraunhofer  region,  the 

concept  of  shadow  boiondary  is  already  meaningless. 

Finally,  well  within  the  geometrical  shadow  region,  0  «9  ,  we  can 

employ  the  expansion  (E.6): 

2 
L(u,v)  »  J^(v)  -  i  ^  J^(v)  -  (J)  J2(v)+  ...     (u  »v),    (8.36) 


which  gives 


^ct 


j^(P9)-i  ^  j^Oe)  -  (^)  JgOe)^ 


(e«-). 

\    r 


(8.37: 


Now  let  us  consider  the  effect  of  the  Fock-type  terms  F(s,t,v)  in  (8.12 ^ 
Just  as  in  (8.21),  we  can  rewrite  (8.1^+)  as 


F(s,t,v; 


Ai( 


Al(x  f       .       -lrt/3    ^  -r  /         -lJt/3 

'-T-f — -r—   exp(-ie   '  sx)J  (v  -  e    ' 

-2irt/3>    -^^  o' 


tx)dx 


xe 


/   ),'■' '  j^      exp(isx)J  (v  +  tx)dx 

J   . .  ,   2in/3  N  o^ 

o   Aifxe   '-^ 


(8.38) 


1/3 
In  the  shadow  region,  for  r  «P  'a  (s  »  1),  it  can  be  shown,  as  in  (8.25), 
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that   isF(s,t,v)  cancels  the  main  term  in  L(u,v),  leaving  as  a  remainder  a  residue 
series  corresponding  to  (5.11). 

For  r  >  p  'a  (s  <  l)  and  0  «p   '   (t  «  l),  the  main  contribution  to 
the  integrals  (8.38)  arises  from  small  values  of  the  argument,  because  of  the 
Airy  functions.   Therefore,  we  may  employ  the  expansion 

(8.39) 

and  similarly  for  the  other  term,  with  the  following  result: 

isF(s,t,v)   =   [f(s)  -  l]j^(v)  -  tg(s)J^(v)  +  ...    (s  <  1,  t  «  1),   (8.I1O) 

where   f (s  )   is   given  by    (8.21)   and 

ilT/3 


J   (v  +  tx)      =     J   (v)    -   txJ,  (v)   + 
o  o  1 


;(s) 


ise 


^  o 


Ai(xe 


Ai(x)  / .       V    ,  -i«/3 

— '  .  '  /-        exp(isx)xdx   -   e 

?lTt/3\ 


;8.1il) 


Ai(xe 


Ai(x)        ,  .  -in/3   \  r, 
— '  .  )^        exp  -le     sxjxax 

-2lrt/3N 


In  particular,  if  also  s  «  1,  we  may  expand  the  exponentials  in  power 
series,  and  we  get,  with  the  help  of  (8.27), 


F(s,t,v: 


^M^  +  iM-j^s  +  ...)J^(v)  -  M^tJ^(v)  + 


(s  «  1,  t  «  1' 


(8.1+2) 


Taking  into  account  (8.11)  and  (8.37),  this  leads  to 
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t. 


.V2 


-; J  exp(ikr  +ik  ^) 

m  0y  2r 


1/3 


2/3 


1  +iM^(|)    f  -M^(|f -'(J)2^  ... 


X  J  (P9)  -  i 


0  2/3^2 
1  +M^(|)    (J)  + 


,1/3 


re 


J^(pe-)  +  . 


(r  »  P  /^a,   G  «  a/r).     (8^3 ! 


Note  that  (8.29)  is  a  particular  case  of  this  result. 

1/3 
In  the  neighbourhood  of  the  axis,  for  r  »  3  '  a,  the  intensity, 

2 
according  to  (8.43),  behaves  like  J  (30)  times  the  intensity  of  the  incident 

wave,  so  that  the  Poisson  spot  actually  corresponds  to  a  "Poisson  cone"  of 

angular  opening 


e  ~  P 


-1 


(8.UU) 


In  this  region,  the  Fock  terms  give  only  a  small  correction,  having  the 

same  angular  dependence  as  the  Fresnel-Lommel  approximation. 

-1/3 
Finally,  if  9  »  P  '  ,  we  must  recover  the  reflected  wave,  given 

by  (6.23).   It  can  be  verified  that  it  arises  from  a  saddle  point  in  one  of 

the  Fock  terms.   However,  we  shall  not  discuss  this  problem  here;  the 

analogous  but  simpler  case  of  the  far  field  will  be  discussed  in  the  next 

Section. 


IX.   The  scattering  amplitude 

A.   Behaviotir  Away  From  Forward  or  Backward  Directions 
For  r  ->  «>,  we  have 


^(r,e)  ^  e^^-"  +  f(k,9)e^^/r. 


(9.1) 
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where  f(k,9)  is  the  scattering  amplitude.   If  9  is  not  too  close  to  0 
or  Ti,    f(k,9)  may  "be  obtained  from  the  expressions  derived  in  Section  VI, 
which  remain  valid  for  r  -^  ooo   Actually,  the  asymptotic  form  (9.1)  is 
already  valid  in  the  "Fraunhofer  region"  defined  hy  r  »  Pa  (cf.  Section 
VIII). 

According  to  (6.19)  and  Fig.  7,  the  angle  w  approaches  zero  as 
r  ^  «j,  so  that  ^  -4  0/2  in  (6„20)  and  the  saddle  point  approaches 

X     =     kp  =  p  cos  I  .  (9.2) 

This  corresponds  to  the  geometrically  reflected  ray  in  the  direction  9, 
as  shown  in  Fig.  11. 

It  follows  from  (6.23)  and  (6.2i|-)  that 


f(k,9)   =  f  (k,0)  +  f   (k,0),  (9.3! 


where 


f  (k,9)   =  -  I  exp(-2ip  sin  f)  (l  + --^-^   +  ...  !     O-K 

^  -   \.    2P  sin-^  I     y 

is  the  "reflection"  amplitude  and 

1^       r.  l/6  -if.     /I    'v — '  exp(iA.  &'  ) 

n    r.\  -irt/6  a,2^    '     ,         ■      a\       J    -iTt/i|  \  n  o^ 

f        (k,9        =     e        '      -(-  rt   sm  9)       <e        '         >  J 

^^^  2  P  ]       .^  [Ai'(-x  )f 

[^  n        n 

00 


+  >  (-1)" 


where,  according  to  (5.7)  and  (5.8), 


exp(iX,  T'  +if)  +  exp(iX,  &'  -  ir' 


r'   =  2mrt  -  9,   6'   =  2m7T  +  9.         (9-6) 
m  m 


*-  z 


Fig.    11 
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The  main  term  of  (9-^)   corresponds  to  the  geometrical  optics 
approximation,  whereas  (9.5)  represents  the  radiation  from  the  surface 
waves. 

The  expression  (9.^)  was  obtained  as  a  limiting  case  of  (6.23), 
so  that  it  represents  the  contribution  from  the  neighbourhood  of  the 
saddle  point  (9.2)  in  (6.6)   or  (6. 31).   It  may  also  be  derived  by 
substituting  directly  in  {6.6)   or  (6.31)  the  expansion 


H^^^(P)  .  (^/   exp 


i  (  P  -  X|  -  ^ 


(9.7) 


which  leads  to 

lAdX.  (9.8) 

r  H-^(p) 


i    r  4^^(p^    (1) 


As  in  (6.6)  or  (6.31),  the  path  P  crosses  the  real  axis  at  the 
saddle  point  (9.2),  at  an  angle  of  -n/h   (Fig.  12).   Since  the  main 
contribution  arises  from  the  neighbourhood  of  this  point,  the  path  of 
integration  may  be  extended  to  infinity  on  both  sides,  provided  that 
the  integral  converges. 

According  to  Appendix  A  and  (C.7),  the  integrand  of  (9.8)  behaves 
at  infinity  like  exp[iA,(2rt-  0)J  in  region  C  of  Fig.  12  and  like  exp(-ixe) 
in  region  A.   Thus,  for  9  /  0,  the  path  r  may  begin  at  infinity  in  C  and 
end  at  infinity  in  A.   In  particular,  it  is  equivalent  to  the  path  r' 
of  Fig»  12,  which  is  taken  to  be  symmetrical  about  the  origin.   This  will 
be  useful  later  on. 

Substituting  (6.11|)  and  (C.7)  in  (9.8)  and  employing  (6.21),  we 


ReX 


Fig.    12 
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are  led  again  to  (9. ^). 

The  result  (9-^)  can  no  longer  be  applied  when  the  correction 
terms  in  the  expression  within  brackets  become  comparable  to  the  main 
term,  i.e.  for  6  <  P    .   Furthermore,  the  asymptotic  expansion  (C.7! 

of  Q^  i(cos  9)  for  0  =  it  -  e  is  only  valid  for  |  A,e  I  »  1  or,  with 

-'""  _i 

X  =   X,  =  P   sin(  e/2)~  p   e/2,    for   e  »  p   2_      ^-^-^^  -^he   conditions  for  the 

validity  of  the   above  results   are 


ri/3 


9  »  P     '^^,  n  -  0     »     p  2.  (9.9) 

In  the  next  sections,  we  shall  examine  what  happens  near  the 
forward  and  backward  directions,  when  these  conditions  are  no  longer 
satisfied. 

B.   The  Neighbourhood  of  the  Forward  Direction 

Let  us  consider  first  the  neighbourhood  of  the  forward  direction, 
defined  by  9  <  P  '  .   In  this  domain,  not  only  does  (9.^)  lose  its 
validity,  but  also  the  residue  series  involving  &'  in  (9- 5)  is  no  longer 
rapidly  convergent,  so  that  we  must  make  a  rearrangement  similar  to  that 
of  Section  VIII,  namely. 


f(k,9)   =  f  (k,9)  +  f   (k,0),  (9.10; 

o        res 

where,  according  to  Section  VIII  and  (C.ll), 


-in/6  B  -"-/^   o  h    v^    m  \~<  exp(2imrtX.  ) 

res'  2'    ^sm  0^    /  ;  X  .  r„  .  .  /    ,-,2  o^  n  " 

-^-r '      ^— '   Ai '  (  -X 


(9.11) 


and  f  (k,0)  is  defined  by 
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o  o 


(r  -^oo),       (9.12) 


where  f     is  given  by  (8.3). 


Since  (8.10)  remains  valid  for  r  ^  "o,  we  find  at  once  that 

-  r       ° 

2    '    2ijt/3    r   Aifxe" 


.    e^^      .fil/3,      9      ,*J   2iK/3    r   Ailxe^ll 


„ . ,    TiTt/s'x 
Ai(xe        '     ) 


X     J 


39  +   (|)        9x1  dx  +  e^'^/^ 


Aifx) 
Ai(xe 


2i«/3  X   '^o 


39  +  (^)     ex 


dx' 


(r  ^c»),  (9.13) 


According  to  (8.3  )j  we  have 


^03  =   ^^^   ^ 


^2°°       P^ 


.(1), 


.X.7t/2 


Hil^(p)P,  i(cos  e)e^^''/  XdX..   (9.14) 


'/ 


/I), 


In  the  integral  from  0  to  (3,  we  can  replace  H^   (p)  by  its  asymptotic 


expansion  for  p  ^  oo^  so  that 


♦o3  *   i  ^  /^x-4(=-  «w^  *  ^t/  ^~"^' 


a^co 


:(1), 


.X.3T/2, 


(9.15) 


X  /    Hi-^^(p)P,  i(cos  9)e^^''/^Xd;.       (r  ^oo). 
Jq  a.      A,-2 


On  the  other  hand,  since 


2inX 


tan(jtX,)   =  i  -  2i 


1  +  e 


2ijtX  ' 


(9.16: 


it  follows  from  (2.3U),  with  a   =  a      (cf.  Fig.  2  and  Fig.  6),  that 
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Og" 


=  (— )  e   '   /  H;  '(p)  P  j^(cos  e)e   '  MX  -  ^ , 


'2p 


(9.17) 


loo 


'^(p,e; 


( — )  e   '   /  h;  '(p)P^  i(cos  e)e   '   


2inX. 


2inX 


Ad  A..  (9.18) 


In  the  last  Integral,  the  path  of  integration  has  been  shifted  to  the 
positive  imaginary  axis,  which  is  allowed  because  of  the  extra  convergence 
factor  exp(2inX,).   Substituting  (9.17)  in  (9-15)^  we  get 


„   .   ikr  p 

ip  cos  0   1  e /  ^   /    ^\,T,    A      /     \ 

t„^  "x  <s                     + /  P,  i(cos  e)>aiX  +  A      (r  ^oo) 

03              k  r  J   X-2 


(9.19) 


The  second  term  of  this  expression  can  be  evaluated  by  inserting 

for  P  j_(cos  0)  the  uniform  asymptotic  expansion  (Cll)  and  integrating 
A,- 2 

term  by  term  with  the  help  of  the  well-known  formula 


-n+1  ^  ,  ,j^        -n+1  ^   ,    . 
X  J  (x jdx  =  -X     J  ^ (x, 

n'  n-1 


(9.20) 


The  result  is 
P 
i 


1. 

-      I    P^  i(cos  0)XdX,  =  ika^(-r^) 
k  J      \-2  sm  0 


m 


(9.21) 


On  the  other  hand,  making  \  =  i|a  in  (9.l8 ),  we  find 


|A|  <   (2^ 


JL   °° 

2   r 


"^"/2|h!1)(p)|  |p.   i(cos  0)1  e-^'^Vn.     (9.22; 
'  IM-         l(-i-2 


According  to  Watson  (2^),  we  may  employ,  for  all  u  >  0  and  p  »  1, 
the  asymptotic  expansion 


„(1),  .     ,2,2  2^  2,-l/l|    1. 
^n-M  (P^  ~  (-)  P  +  1^  )   '  exp  ^1 


'p  +  |i  )^  -  n  sinh'  -^  -  i 


P  "  I 


.f 


>  , 


J 

(9.23) 


so  that 


Furthermore,  according  to  (C.12), 


Substituting  (9.2^)-)  and  (9.25)  in  (9-22),  we  find 

(r  ->  oo). 


*   _-l 


|A|     <     -i^   (cos   0)   2   p 
9it 


■1/3 


(9.21^) 


p.      1  (cos   0)1      <     (cos   0)"^  e"^*^^  (O<0<rt/2).        (9.25) 

11-1-2 


(9.26) 


Since  we  are   only  interested  in  the  domain  0  <  3      '    ,    it   follows 
from   (9.19),    (9.21)  and   (9.26)   that 

ikr 


i 


03 


,ip  cos  e  ,  ^^^^i  rhj^ 

^sm  0     30 


+  e^(3 


(r  ^-).   (9.27; 


Combining  this  with  (9.13  )j  we  finally  obtain 

•1  2   „   1 
ika   ,  0   \p 


fo(k,0) 


(±-^f   [g^,(k,0).g^(k,0)],    (9.28) 


where 


^^(k,©)   =  2J^(P0)/p0, 
g^(k,0)   =  {2/iif^^F{0,t,Y), 


(9.29) 


(9-30) 


F(s,t,v)  being  defined  in  (8.IU). 

Except  for  the  substitution  sin  0^0,  which  is  allowed  in  this 
order  of  approximation,  (9.29)  coincides  with  the  classical  Fraunhofer 
diffraction  pattern  of  a  circular  disc  or  aperture  (26).   The  other  term 
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g  may  be  called  the  Fock  correction  term  (27). 
F 

The  forward  scattering  amplitude  is  obtained  by  setting  9=0 

in  (9.28): 

f^(k,0)  =  |ika^(l  +  2Cp'^/3+  ...),        (9.31) 

where  C  is  given  by  (7.^2)  and  (7.^3).   This  result  also  follows  from 
(8.29),  by  taking  z  »  |3a. 

The  total  cross  section  is  obtained  from  (9.31)  with  the  help  of 
the  well-known  "optical  theorem" 

a  =  ^  'S^f{^,0)     =  rta^(2  +  1.9923P"^^^  +  ...),   (9-32) 

where  we  have  employed  (7.^3)* 

This  result  was  first  derived  by  Rubinow  and  Wu  (20).   Higher-order 

-2/3 
terms  in  the  expansion  in  powers  of  P  '   have  been  computed  by  Wu  (2^) 

and  by  Beckmann  and  Franz  (28).   They  involve  the  coefficients  M  defined 


in  (8.27). 

If  0  «  P  '    ,   we  may  employ  {Q.h2.)   with  s  =  0.   The  resuli 

1/3 


h  \    J,(pe)   ,  2/3 


f^(k,0)   =  iika2(-f-^)  U   -Vi-  -  (f ) 


M^J^(P0)  -  M^(|)   0J^(P9) 


=  WJa^^-^ 


2     )     P0 


•+  e^^/^p'^/^ [1.9923  Jo(Pe; 


+  0.6706  pl/39  j^(P9)  +...]  I   (0  «  P"^/3)^      (5^33) 


where  we  have  employed  (8.28).   This  is  also  a  limiting  form  of  (8.i<-3)  for 
r  »  Pa. 
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The  first  terra  of  (9.33)^  which  gives  rise  to  the  well-known  forward 

-lA 
diffraction  peak,  dominates  throughout  the  region  0  «  p   '  .   The  corrections 

are  of  the  order  3  '  0. 

-1/3 
If  0  becomes  »  P  '  ,  we  must  recover  (9- 3).   To  show  this,  let  us 

rewrite  F(0,t,v)  as  follows  (cf.  (8.lU)  and  (8.21)): 

00 

-iTr/3  _,._    ,       s  r  Aifx)  -,   ,  .    -In/s    \-, 

F(o,t,v)    =    /   ---z^T^r   Jo^^  -  *^       ^^<^^ 

^   Ai ( xe    '  ) 
o    ^ 

(9.34) 

+  r -h^H:r     J  (v  +  tx)dx  . 

J      Ai(xe2W3)   o^ 
o 

Let  us  make  the  decomposition:   J  =  [H    +  H   ]/2  and  substitute 

■^  o  o  o       ' 

X  =  e  '    x'    in  the  first   integral  in  H        ,    employing  also  the  relation   (D.3). 

The  result   is 

00 

2e-i«/3i.(0,t,v)   =    r ^^^7-  H^2)(,   .  te-'^/3x)dx 

^      Ai(xe  '-') 


"j^T.     H(2)(,.tx)dx./-;^iifi73^H(^'(v.tx).x     (9.35) 

2iTt/3 


^      Ai(xe        '     )  ^     Ai(xe 

o  r         ^ 


-in/3 


-^  V-  /         ^^^\^)^^ , 


V 


where  P  is  the  path  of  integration  shown  in  Fig.  10.   This  path  can  be  closed 
at  infinity,  reducing  the  integral  to  a  residue  series 


^  ' — '  exp(ie    tx 


J  \W        Z_i   [Ai'-(-x  )f 


(9.36) 


(-X  )] 
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where  -x  are  the  zeros  of  the  Airy  function  and  we  have  replaced  H    by  its 
n  o 

asymptotic  expansion. 

Since  v  »  t  »  1,  the  Hankel  functions  may  be  replaced  by  their 

asymptotic  expansions  also  in  the  other  terms  of  (9-35).   By  partial  integration, 

we  find  / 

00  e   '  i 

f  H^l)(x)dx  =  e^'^/^d;)'   e^M  1  -.  ^  (v"!)]  .       (9.37) 

V 

The  asymptotic  expansion  of  the  second  term  of  (9.35)  can  also  be 
obtained  by  partial  integration.   The  first  term  has  a  saddle  point  at 

>/3  +2, 


X  =  e 


tyk,  (9.38) 


as  can  be  verified  by  replacing  the  Airy  functions  by  their  asymptotic  expansion 
(D.  ij-).   In  addition  to  the  saddle  point,  we  must  take  into  accoimt  the  contribution 
from  the  lower  limit  of  integration.   The  result  is 


00 


O 


Ai  X  „(2  ,    ,  -irt/3  ^^ 

Xt-'-j--  H^  (v  -  te   '  X  dx 

,.,      -2lJt/3^   o  ^ 
Ai(xe    '  ) 


Ai(xe'"^"'-^) 


^1^  Hf'(v.t.,.x»  (f )*e-/%  (.,.  .  li^) 


o 


2  \^    e 


\   .-i«/3 


^(  ^j  — T"  ^^P 


(9.39) 


where  the  first  term  represents  the  contribution  from  the  saddle  point  and  the 
second  one  those  from  the  lower  limits  of  integration. 

It  follows  from  (9.35)  to  (9.39)  that 
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F(0,t,v) 


tN* 


2v 


■) 


exp   ( -iv  + 


it!" 

12 


1  ^  2   ^2 


cos      V 


ijr/l2     r— 1     exp(iv  +  ie     '    tx    " 
e     '  \  n 


(2jtv; 


^— — 1  exp (IV 

^   [Al'( 
n 


(v  »  t  »  i; 


(9.to) 


[Al'(-x^)]' 


Substituting  this  result  in  (9-28)  to  (9.30),  and  replacing  J  (P0)  by 
its  asymptotic  expansion,  it  is  found  that  the  main  term  of  this  expansion  is 
cancelled  by  the  second  term  of  (9.^0),  leaving  us  with 

3n^ 


fo(k,e) 


a  ( 0. 

2  Vsin  0 


ejrp  <-2iP 


3:  V2 


{9.hl) 


+  e 


-5ijr/l2 


\i/« L 


{n   sin  0)2 


exp(iA,  0) 


n 


[Ai'(-x  ) 
^   n 


s— -  (0»p-i/3; 


The  second  term  is  identical  to  the  residue  series  involving  &'  m  (9. 5). 

o 

The  exponent  of  the  first  term  is  the  expansion  of  -2ip  sin  p  up  to  the  term 

3  3 

of  order  |30-^.   Thus,  in  the  domain  where  P0  »  1  but  the  next  term  in  this 

expansion  can  be  neglected  (and  at  the  same  time  0/sin  Q  v  1) ,   the  first  term 

of  (9.^0)  coincides  with  the  main  term  of  the  reflection  amplitude  (9. ^). 

-1/3 
Therefore,  (9. 10)  goes  over  smoothly  into  (9.3)  when  0  becomes  »  3 

but,  just  as  was  found  in  Section  VTI,  the  Fock  functions  cannot  be  employed 

for  too  large  values  of  9.   Their  angular  domain  of  validity  is  just  sufficient 

to  make  a  smooth  transition. 


The  transition  from  the  forward  diffraction  peak  to  the  region  of  geometrical 

-1/3 
reflection  takes  place  in  the  domain  0  ~  P  '  .   In  this  domain,  we  must  employ 

the  expression  (cf.  (9.10),  (9-11),  (9.28)  to  (9.30)  and  (9.3^)): 


-  74 


f(k,e)        2  ik=   (^si„  ej    -^2     pe     +1^,3;      [J  2i;73 


o  V^"  -  "  '"^'  (§)'  '»;  -  -  J    1,7^.-/3, 


X    J    (pe  -  e""^'/^  (^)      ex)  dx  +     '    — ^^^' 


if  ax)  ax].... I    .e-W6,(|)^/' 


X    J^  (  3G  + 


]^  CO  ,  . 

m=l  n  ^      n 

(O<0<p"^/3)^  (9.te) 

-1/3 
where^  for  0  ~  P  '    ,   the  Fock-type  functions  should  be  computed  "by  numerical 

methods.   Some  related  ftinctions  have  already  been  tabiilated  (22,29),  tut  there 

seem  to  exist  no  tables  for  those  appearing  in  {9.k2), 


C.   The  Neighbourhood  of  the  Backward  Direction 

There  remains  to  consider  only  the  neighbourhood  of  the  backward  direction, 

i.e.,  according  to  (9.9),  the  domain 

1 
©  =  It  -  e,  e  <  p'-  .  (9.^3) 

The  expression  for  f(k,9)  in  this  domain  may  be  obtained  similarly  to  (9.8), 
by  substituting  (9.7)  in  (6. 30)  and  (6.35): 

f(k,0)  =  f^(k,9)  +  f^^^(k,0),  (9.U) 

where 
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f  (k,e)   =  ^   /    -frr P  i(-cos  e)tan(jrA.)e  ^"^Vx.,  {9.h^] 


A/ 


Ex  r  O) 
^^-^— -  P,   (-cos  9)  .  (9.46: 

cos  nX         X     I  \^        . 


n 


The  path  of  integration  in  (9.^5)  is  the  upper  half  of  the  path  P'  shown 
in  Fig,  12.  Actually,  (9.^5)  can  also  be  obtained  from  (9.8),  by  taking  the 
integral  along  the  symmetric  path  P'  and  then  applying  the  same  transformation 
that  led  from  (6.30)  to  (6. 31). 

Substituting  (2.18),  (3.16''.  and  (9.^3)  in  [9.^6),    and  employing  the 
uniform  asymptotic  expansion  (C.ll),  we  get 

1/3  '    ^   — 


res-  \2J         \3in  c/   ^/_,       /^j 


m=0       n 


exp[i(2m  +l)rt?i.  ] 


{9.hl\ 


[Ai'(-xJ 


2     o'  n 


which  is  also  a  limiting  case  of  (6.38).  In  particular,  if  e  »  P  ,  we  may 
replace  J  ( X  e )  by  its  asymptotic  expansion  and  (9.^7)  goes  over  into  (9-5)^ 
as  may  easily  be  verified." 

In  (9.ij-5),  we  may  employ  the  expansion  (6.1^,  which  remains  valid 
even  for  \x\    «   3.   Acoording  to  f9.2),  the  main  contribution  to  the  integral 
arises  from  the  neighbourhood  of  the  lower  limit,  i.e.  from  the  domain 

\X\     <    p"*  .  {9.h8) 

Thus,  we  can  expand  the  exponent  and  the  other  terms  of  {6.1h)   in  powers  of 


t6 


A,/P,  with  the  following  result 

e-^''h[^\^)M^\^)      =     ie-2iPexp(-ir/P)  [l  +  -^  -  ^  +  . .  .J   ,      {9.h9) 

where  we  have  kept  all  correction  terms  up  to  the  order  P   ,  according  to  (9.  ^). 
On  the  other  hand,  P  j^(-cos  0)  =  P  ^.(cos  e)  and,  according  to  (9-^3) 

A,"  2  ^"  2 

and  {S.hQ),   the  relevant  portion  of  the  domain  of  integration  corresponds  to 
I  \e  [  <  1,  so  that  we  may  employ  the  expansion  (C.9)« 

Making  the  change  of  variable  appropriate  to  the  steepest  descent  path 
(cf.  Fig.  12) 

X     =     e^^^'/^p^x   =ax,  (9.50) 

we   finally  get   from    (9.^5) 


f   (k,«-  e)      =     iae'^^*^ 
r 


where 


[  1  +  -rjr  )      /   exp(-x    )J   (wx)tan(jtax)xdx 
o 

00 

i         r  ,      2.      ,      s        ,  \5,  uj.2e 

+     r-pR      /      exp(-x    )J   ((Jx)tan(3ta  x)x'^dx  +  r  sm     - 


r  2  2  2    e 

X      /   exp(-x    )J   (wx)tan(nax)x  dx  -   sin     -  (9.51) 


/2  1  2   £ 

exp(-x    )Jp(wx)tan(jt  ax)xdx  +^    sin     - 


00 

X      /   exp(-x    )J^(wx)  tan(:tax)dx   +  Q' {?>'    )      , 


u)     =     2e^^"'^  p2   sin  |     =     2a  sin  |  (9.52) 
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and  the  upper  limit  of  integration  has  been  extended  to  <»,  since  large  values 
of  X  give  no  significant  contribution.   Note  that  |w|  <  1  according  to  (9.^3)« 

The  evaluation  of  the  integrals  appearing  in  (9.51)  is  taken  up  in 

2  2  e 

Appendix  F.   The  result  is  given  by  (F.13).   Expanding  (J  =  -4iP  sin  p  in 

2 
powers  of  e  and  taking  into  account  (9.^3)j  we  finally  get 


f^(k,n-  e)  =  -  2  exp 


-2iP  fl  -  4 


(0  <  e  <  p  2)^     (9.53) 


As  may  readily  be  verified,  this  result  coincides  vith  the  expansion  of 

2  -- 

(9.^)  in  powers  of  £  ,  within  the  domain  e  5  P  ^.   Thus,  (9>^)  is  uniformly 

valid  up  to  0  =  n.      This  had  often  been  assumed  in  previous  work,  on  account  of 

the  regularity  of  (9.^)  up  to  9  =  «.   Obviously,  however,  such  regularity 

constitutes  no  proof  of  the  validity  of  (9.^),  since  the  asymptotic  expansion 

(C.7)^  employed  in  its  derivation,  is  no  longer  applicable  in  this  domain. 

Since  (9.^7)  also  remains  valid  for  e  »  P   ,  we  conclude  that 

f(k,e)      =     -  I  exp(-2ip   sin  |)  (  1   + 


2      ^^  2'V  23   3in3| 


/  1  °° 

-^  m=0  n 


exp[i(2m  +1)-kX  ] 

X ^  jjx  (« -  0): 

[Ai'  (-x^)] 


uniformly  throughout  the  whole  domain 

■1/3 


p-^/-5     «    0      <    Tt  .  (9.55) 
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Together  vith  (9.h2),   this  determines  the  behaviour  of  f(k,0)  for  0  <  0  <  it. 

D.   Direct  Transformation  of  the  Scattering  Amplitude 

The  expressions  for  f(k,0)  employed  in  the  above  discussion  were  obtained 
by  letting  r  -^  <»  in  the  representations  previously  derived  for  t(r,0).   The 
advantage  of  this  method  is  to  make  clear  the  connection  between  the  behaviour 
of  the  wave  function  in  the  near  and  in  the  far  regions,  as  well  as  the  physical 
interpretation  of  the  various  terms.   However,  one  may  ask  whether  it  is 
possible  to  bypass  the  limiting  procedure  and  to  derive  the  same  resiolts  directly 
from  the  partial-wave  expansion  of  f(k,0), 

00 

f(k,9)   =  ^^^^    [S^(k)  -  l]P/cos  0),  (9.56) 

i=0 

where  S-  is  given  by  (2.3). 

It  is  well  known  (30)  that  Watson's  transformation,  in  the  form  usually 
applied  to  Yiikawa -tjrpe  potentials,  cannot  be  applied  to  (9.56)  in  the  case  of 
a  cut-off  potential  (including,  in  particular,  the  case  of  a  hard  sphere).   This 
is  essentially  due  to  the  asymptotic  behaviour  of  the  S-function  as  [  X,|  ->  °°, 
(arg  X,|  -^  n/2. 

Modified  versions  of  Watson's  transformation  have  been  employed  for 

this  purpose  (31>  32).   However,  the  proposed  modifications  still  lead  to 

5 
singular  integrals  ,  so  that  the  diffic\ilty  is  not  overcome. 

It  will  now  be  shown  that  it  is  indeed  possible  to  derive  all  the 
representations  for  f(k,0)  employed  above  directly  from  (9.56). 
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Let  us  start  by  applying  Poisson's  sum  formula  (2.12): 

00  oo 

f(k,G)  =  I   2  (-1)""  J    ^1  -  S(X.,k)]  P^.i(cos  e)e2^"^''V\,    (9.5?: 
in=-oo       o 

where  S(X,k)  is  given  by  (3.1).   It  follows  from  (2.15)  that 


S(-A.,k)   =  e^^"^  S(X,k),  (9.58) 


so  that,   making   X,  ^  -X  in  the   sum  from  m  =  -1  to   -<»,   we  get 

o 
f(k,e)      =     i     2   (-1)"^    i     f    pi'^^  -   S(X,k)]p^_i(cos   9)e2i'"«Vx 


m=0 


[l   -   S(X,k)]P^_i(cos   0)e^^"''Vx^     .  (9.59; 


Note  that 


.^i'^^  -    S(X,k)      =     2e^"^J_,(P)/H[^\p).  (9.60) 


According  to  Appendix  A  and  (C.8),  the  integrand  of  the  first  integral 
in  (9.59),  for  all  m  >  0,  goes  to  zero  at  least  exponentially  for  |  \|  -^  «> 
in  the  second  quadrant.   Thus,  we  may  shift  the  path  of  integration  to  the 
positive  imaginary  axis,  from  ioo  to  0.   Writing 

e^^''^  -  S(X,k)   =  e^^"^  -  1  +  1  -  S(X,k) 

in  the  corresponding  terms  for  m  >  1,  and  employing  the  identity 

J  (-l)"^Ce2i"^  -  1)  e^i"^'^^  =  e^^"^  ^   ^  ^.^^  -  ^  ,   (9.6l) 
m=l  ^  '^ 
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we  get 

f(k,G)    =    f^(k,e)  +  f^g3(k,e),  (9.62) 

where 

00 

f   (k,G)      =7      /   [1   -   S(A.,k)]P.    i(cos   9)}^X 

O  ^      J  X~2. 

O 
.  /? 

-|      /  S(X,k)P^_i(cos   9):\d\  +  A^  ,  (9.63) 

i  00 

„.       8  2irtX, 

A.      =     ^      /  ^   ^.        P^    i(cos   Q)MX,  (9.64) 

1  k     ,/  ,         2irtX     X-2  ' 

y  1  +  e 
100 

00 
f^^g(k,0)      =     I     ^  ^"^^""Z    ^^   "   S(X,k)]P^_i(cos   G)e2^"'''^XdX,       (9.65) 
m=l  C 


and  C  is  the  path  shown  in  Fig.  13,  going  from  ioo  to  0  and  from  0  to  00  . 

It  follows  from  Appendix  A  that  this  path  may  be  closed  at  infinity 

in  the  first  quadrant,  so  that 

00         00 

f   (k,G)   =  T^   >   (-I)'"   >   X.  r  exp(2imrt\  )  P,   i(cos  G),    (9.66) 
res^  '    '  k  /  \  ^  /   I  n  n  -^^     n   X.  -2 

m=l       n=l 
where  X,  are  the  poles  of  S(A,,k)  and  r  are  the  corresponding  residues,  given 
by  (3.1'+).   This  corresponds  exactly  to  (9.11). 

On  the  other  hand,  we  may  rewrite  (9.63)  as 

f,(k,G)  =  f^^  +  fQ2  ^  ^03'  ^9.67) 

where  >. 

^01  ^  ^02  =  J  (  /   -   /  )  ^^£|^   M^-^  ®)^^ 

\i»     o  J   h[i)0) 

^      00 

+  2i   r   yP)   p   ^(cos  G)AdX  (9.68) 


ImX 


ReX 


Fig.    13 


and 
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P 

f^^  =  7   /  P,  i(cos  e)XdX  +  A  .  (9.69) 

03     k   J    X-2  J- 


The  first  integral  in  (9.68)  may  also  be  taken  along  the  path  C'  from 
a  00  to  P  shown  in  Fig.  13 .   Since  the  main  contribution  to  both  integrals  arises 
from  the  neighbourhood  of  X  =  P,  where  we  may  employ  the  approximations  (7»33) 
and  (7.3^1),  we  see  that  (9-68)  corresponds  to  (9.13).   Similarly,  (9.69) 
corresponds  to  (9.19).   Just  as  in  (9.26),  it  can  be  shown  that 

1^.  I  <  — o— ^ T  ,  (9.70) 

9rt  k(cos  e)2 

so  that  the  contribution  from  ^  ,  which  is  independent  of  a,   would  be  included 
in  the  correction  term  of  order  P   in  (9.27). 

Thus,  (9.62)  leads  to  the  same  results  as  (9.10)  and  is  the  appropriate 
splitting  of  f(k,e)  in  the  domain  9  <  P  '^. 

Now  let  us  consider  the  domain  0  »  P  '  .   In  this  case,  the  integrals 

containing  Q^  2(cos  9)  in  (9.63)  can  also  be  reduced  to  residue  series.   For 
A.- 2 

this  purpose,  let  us  add  and  subtract  the  (convergent)  integral 

o 
/  Q^?i(cos  e)MX 

loo 

and  rewrite  (9-62)  to  (9.65)  as  follows: 

f(k,9)   =  f^(k,9)  +  f^^Jk,9),  (9.71) 

where 

00 

f   (k,9)  =  f   (k,0)  +  (2jt/k)  >   X.  r  qJ^-*!  (cos  9),         (9.72) 

res^  '        res^  ■*  '    \   /  '  / ^  n  X,  -2 

^ — t'  n 

n=l 
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f    being  given  by  (9.66),  and 


res 


[^hcos  e) 

A,-2 


f^(k,e)  =     -  ^    J    rs(\,k)Q^^i.(cos  9)  +  Q, 
ioo 

00 

+  ^  /  [l  -  S(A.,k)]Q^^i(cos  e)AdX+A^. 


MX 


(9.73) 


Prom  the  discussion  given  in  connection  with  (9.8)  and  (9.59)  and 
from  (C.7)  it  follows  that  the  path  of  integration  in  the  first  integral 
of  (9-73)  may  be  deformed  away  from  the  imaginary  axis  into  the  second 
quadrant,  so  as  to  coincide  with  the  upper  half  of  the  path  P'  in  Fig.  12. 
Similarly,  according  to  Appendix  A  and  (C.7),  the  path  of  integration  in  the 
second  integral  may  be  deformed  within  region  A  of  Fig.  12,  so  as  to  coincide 
with  the  lower  half  of  P' .   This  leads  to 

-000  o 

f^(k,e)  =  -i  r  s{\,k)q,^^_i{cos Q)Mx - ^  r 


O  oo 


O  00 


O.lh) 


q,[^licos  e)  +  Q^^\(cos  e)-  2  ^  p.  ,  .,  ^ 

1_  A.-2  ~A.-2  2.+ e 


P  j_(cos  0) 


MX, 


where  the  integrals  are  taken  along  the  path  P' .   We  have  also  made  A,  ->  -X.  in 
the  integral  of  Q,\    I     and  employed  (9.64)  (with  the  path  of  integration  shifted 

A- 2 

to  P'  ). 

It  follows  from  (C.5)  and  (C.6)  that  the  expression  within  square 
brackets  in  (9.7^)  is  identically  zero,  so  that  we  are  left  with 


— O  oo 


f^(k,e)  =  -J    J  S(\,k)Q^^](cos  0)XdX, 


(9.75) 


(J  00 


in  exact  agreement  with  (9.8).   Since  (9.72)  also  corresponds  to  (9-5),  we 
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see  that  the  splitting  (9- 71)  is  equivalent  to  (9.3  )• 

Furthermore,  just  as  in  (9.^5)^  we  may  rewrite  (9- 75)  as 

a  00 

f  (k,9)   =  -  7   /  S(X.,k)P.  i(-cos  e)tan(nA.)e"^'^^AdA.      (9-76) 
r  ^      J  ^~2 

o 

and,  with  the  help  of  the  identity  (6.33)^  we  may  rewrite  (9-72)  as 

00  00 

f^^Jk,e)   =  -(2rti/k)  2_j  (-1)™  2  Vn^^I^^^^  ""  ^^""^J^^  _^(-cos  e), 

m=0       n=l  ^ 

{9.11) 

which,  according  to  (2.l8),  is  equivalent  to  (9- ^6). 

The  results  (9-76)  and  (9.77)  remain  valid  up  to  0  =  rt  and  correspond 

to  those  employed  in  (9.^^).   In  this  form,  the  splitting  (9.71)  may  be 

-1/3 
employed  for  P  '-^  «  0  <  n. 

In  conclusion,  according  to  (9.62),  (9. 61|) ,  (9. 66),  (9.67)  to  (9.69), 
(9.71),  (9.76)  and  (9.77),  the  modified  Watson  transformation  may  be  expressed 
as  follows: 


f(k,e)      =     -  J      /     S(A.,k)P^_i(cos   9)Xd\  +  J     /    [l   "   S(X.,k)] 
C  P 

P  o 


„.     ,.  2irtX, 

2i    /  e 


X  P^_i(cos   9)XdX  +  I    J     P^_i(cos   e)XdX.  +  fj     -^-^~l  (9.78) 

00  00 

X  P     i(cos   0)AdX  +  ^      >      (-I)'"    >      ^  r  exp(2imrt\   )P,      i(cos   ©), 


m=l  n=l 

where  the  path  C    is   shown  in  Fig.    I3,   and 


-  Qk 


f(k,e)   =  r   /  S(A,,k)P  i(-cos  e)tan(jrA.)e"^"^AdA. 


^"  (9.79) 

OO  00 

m=0       n=l 


where  the  integral  may  also  be  expressed  in  the  equivalent  form  (9.75). 

Both  of  these  representations  are  exact  and  their  terms  have  a  direct 
physical  interpretation,  as  discussed  above.   They  may  be  employed  to  obtain 
higher-order  corrections  to  the  results  derived  in  the  present  section: 
(9.78)  should  be  employed  f or  0  <  0  ~  p""""'-^  and  (9-79)  for  p,~-^'^  «  Q     <  n. 


X.  Conclusion 

The  main  results  obtained  in  the  present  work  may  be  siamnarized  as  follows; 

(i)  The  high-frequency  behaviour  of  the  wave  function  in  scattering 
by  a  totally  reflecting  sphere  may  be  completely  determined,  both  in  the  near 
and  in  the  far  regions  of  space,  by  means  of  a  modified  Watson  transformation, 
based  upon  Poisson's  sum  formula.   Each  term  in  the  transformed  series  has  a 
direct  physical  interpretation.   This  procedure  has  the  advantage  that  it  does 
not  require  a  reevaluation  of  the  whole  residue  series  in  going  over  from  the 
shadow  to  the  lit  region,  but  only  of  that  part  of  the  lowest-order  term  that 
would  not  correspond  to  a  rapidly  converging  series,  thus  showing  clearly  the 
connection  between  shadow  and  lit  region. 

(ii)  It  is  necessary  to  apply  different  representations  in  the  forward  and 

backward  half -spaces.   It  is  already  clear  from  the  singularity  of  P.  i.(x)  at 

A,- 2 
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X  =  -1  that  one  cannot  have  a  single  representation  for  all  values  of  0:   one 

needs  a  representation  in  terms  of  P,  j^fcos  0)  near  0=0  and  one  in  terms  of 

A.- 2 

P,  i.(-cos  0)  near  0  =  it.   The  present  treatment  is  based  upon  the  integral 
A,- 2 

representations  (2.3^)  and  (2.35)  of  the  primary  wave.   In  the  forward  half -space, 
we  employ  (2.l6)  and  (2.17);  in  the  backward  half-space,  we  must  employ  (6.29), 
(6.30)  and  (6.3^). 

(iii)  A  rigorous  proof  of  Watson's  transformation  is  based  upon  a  study 
of  the  behaviour  of  the  integrand  for  •[  x\    -^  oo,  arg  X  -^  5t/2,  as  well  as  of  the 
asymptotic  behaviour  of  the  poles  and  their  residues.   It  follows  from  this 
discussion  that  the  residue  series  (k.l'^)   converges  in  the  whole  forward 
half -space  0  <  ©  <  Jt/2,  but  it  is  only  useful  in  this  form  in  the  domain  where 
it  is  rapidly  convergent,  i.e.  within  the  shadow  region.   This  happens  in  the 
domain  r«|3'a,  ©  -0»P   ',  where  0  =  sin   (a/r)  is  the  shadow  boimdary 
angle.   This  domain,  which  may  be  called  the  deep  shadow  region,  is  represented 
by  the  shaded  area  in  Fig.  l4. 

(iv)  The  wave  function  in  the  deep  shadow  region  is  given  by  (5«7), 
which  represents  a  superposition  of  "diffracted  rays"  arising  from  the  siorface 
waves  associated  with  the  poles  of  the  S-function  in  the  complex  angular 
momentum  plane.   These  poles  do  not  show  the  typical  Regge  behaviour  associated 
with  Yukawa -type  potentials.   The  physical  interpretation  of  these  terms  is 
in  agreement  with  Keller's  geometrical  theory  of  diffraction:  they  give  rise 
to  an  exponentially  damped  wave  f "unction  in  the  angular  variable.   However, 
this  interpretation  can  be  applied,  at  a  given  frequency,  only  to  the  lowest-order 
poles,  and  loses  its  validity  when  the  damping  within  one  wavelength  becomes 
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appreciable.   The  higher-order  poles,  however,  give  a  negligible  contribution. 
The  exponential  damping  in  the  deep  shadow  shoiild  be  contrasted  with  the  much 
weaker  damping  found  in  the  case  of  a  circular  disc:  the  shadow  of  the 
sphere  is  much  darker  than  that  of  the  disc  (27,  p.  ^63).   Physically,  the 
reason  for  this  result  is  that  the  intensity  thrown  into  the  shadow  by 
diffraction  at  the  edge  of  the  disc  is  much  greater  than  that  due  to  diffraction 
aroimd  a  curved  surface.   Classical  diffraction  theory,  which  predicts  the  same 
behaviour  for  the  sphere  and  the  disc,  fails  in  this  domain. 

(v)   In  the  lit  region,  sufficiently  far  from  the  shadow  boundary 

-1/3 
(0-0  »  P    ),   the  WKB  expansion  for  the  wave  function  has  been  confirmed 

up  to  the  second  order.   The  expression  (6.23)  '£ot   the  reflected  wave  is 

valid  both  in  the  forward  and  in  the  backward  half -space,  although  the 

derivation  is  different  in  the  two  cases.   In  this  region,  we  also  find  the 

continuation  of  the  surface  waves,  but  they  are  masked  by  the  much  greater 

contribution  from  the  incident  and  reflected  waves,  except  in  the  immediate 

vicinity  of  the  shadow  boundary. 

(vi )   On  the  surface  of  the  sphere,  Kirchhoff's  approximation  (7»l)-(7.2) 
for  the  normal  derivative  of  the  wave  function  is  accurate,  except  within  the 
penumbra  region,  |0  -  -|  <  P  ' -^    (Fig.  1^).  The  behaviour  in  this  region  is 
described  by  Fock's  function  (7. 11  )-(7. 13  ) j  which  interpolates  smoothly  between 
the  values  on  the  lit  and  on  the  shadow  regions.   However,  it  cannot  be  employed 

too  far  beyond  the  penumbra.   The  penumbra  or,  equivalently,  the  corresponding 

-1/3 
angular  momentum  domain  | A,  -  P |  <  P    ,  is  responsible,  in  the  sense  of 

Huygens'  Principle,  for  the  main  corrections  to  classical  diffraction  theory. 
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(vii)  The  neighbourhood  of  the  shadow  boimdary,  |0  -0  [  «  p  '  ,  at 
distances  P  'a  «  z«P   a,  is  denoted  as  the  Fresnel  region  in  Fig.  14. 
In  this  region,  the  main  term  of  the  angiilar  diffraction  pattern  corresponds 
to  the  classical  Fresnel  pattern  of  a  straight  edge,  with  small  corrections, 
representing  the  effect  of  the  cuirvatTare  of  the  sphere  (cf.  (7.^^)).   One  of 
these  corrections  is  the  shift  of  the  shadow  hotindary,  denoted  by  s  in  Fig. 
lU.   It  is  given  by  (7.50),  in  agreement  with  a  conjectiAre  made  by  Rubinow 
and  Keller  (21), 

The  solution  in  this  region  is  still  not  in  agreement  with  the  result 
given  by  classical  diffraction  theory,  which  would  be  analogous  to  the  Fresnel 
pattern  of  a  slit,  rather  than  an  edge.   In  fact,  it  would  contain,  in  addition 


(2) 
to  ^   ,  which  is  actually  replaced  by  a  rapidly  convergent  residue  series 


to  (7.28),  a  contribution  from  the  diametrically  opposite  edge,  corresponding 

(2) 

to  ^^  '  ,   wh 

^o  ' 

(cf.  (7.21)). 

(viii )  In  the  Fresnel -Lommel  region,  0  <  0  ,  P  '  a  «  r  «  Pa  (Fig.  1^), 
the  main  term  of  the  wave  function  corresponds  to  Lommel 's  classical  solution 
(8.17)  for  the  diffraction  of  a  plane  wave  by  a  circular  disc.   The  main 
correction  term  is  given  by  the  Fock-type  function  (8.IU),  which  also  gives 

rise  to  a  smooth  transition  to  the  deep  shadow  and  to  the  lit  region.  Along 

1/3 
the  axis,  starting  at  a  distance  z  ~  P  'a,  we  find  the  well-known  Poisson  spot, 

which  actually  corresponds  to  a  cone  (Fig.  Ik)   of  angular  opening  ~P   , 

surrounded  by  diffraction  rings  (cf.  (8.^+3)). 

(ix)   In  the  Fraunhofer  region,  r  »  Pa,  the  wave  function  is  given  by 

(9.1).   For  p'  '^«   0  ^  «,  the  scattering  amplitude  is  given  by  (9-5^),  the 

-1/3 
main  term  of  which  corresponds  to  geometrical  reflection.   For  0  «  p  '  ,  the 

amplitude  is  dominated  by  the  forward  diffraction  peak  (cf.  (9.33 ))>  which 


corresponds  to  the  classical  resiilt  for  a  circular  disc.   In  the  transition 
domain  0  ~  P  '  ,  we  have  to  employ  (9.^)^  which  again  may  be  continued 
smoothly  up  to  the  region  of  geometrical  reflection.   Higher-order  corrections 
may  be  obtained  from  the  exact  representations  (9.78)  and  (9.79)>  which  result 
from  applying  the  modified  Watson  transformation  directly  to  the  scattering 
amplitude. 

(x)  Fock-type  functions  such  as  (7-13),  (7-39),  (8.1U),  (8.20)  and 
(9-3^)  play  an  important  role  in  linliing  the  domains  of  geometrical  optics 
and  classical  diffraction  theory.  In  view  of  this  role,  which  probably  is 
not  restricted  to  the  present  example,  but  is  of  more  general  validity,  it 
would  be  desirable  to  construct  tables  and  graphs  of  these  functions.  Only 
in  a  few  cases  is  this  material  presently  available. 

Possible  applications  and  extensions  of  the  present  treatment  include 
the  problem  of  a  transparent  sphere  (square  well  potential  in  quantum  mechanics), 
which  will  be  discussed  in  a  subsequent  paper. 
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Appendix  A.   Asymptotic  behavioiir  of  the  cylindrical  functions 
The  asymptotic  hehavioTir  of  the  cylindrical  functions  Z  (x),  x  >  0, 

A. 

in  the  complex  X   plane  may  be  derived  from  the  formulae  given  by  Watson 
(25,  p.  262).   The  results  are  graphically  presented  in  Fig.  15-   The 
notation  is  as  follows: 

A(X,x)      =      (2AnA.2  -  x2)-lA,  (A.i; 


2    2 


a{X,x)      =      (X     -  X    r    -   XjUv 


A  +  i^ 


2  2x3 


X        X 


^], 


(A.  2; 


where  the  branch  of  ( X,^  -  x^)^  to  be  taken  is  specified  by  the  condition 


.2    2A 


{->?   -   x^)2  -^     X     =      \X\    exp(i?)  (-rr  <  cp  <  n)        for  \x\    ^  «>  .   (A. 3) 


Thus, 

1 


A  ^(ttT.     e"-  ^  (If^      for    |x|^-.  (A.U) 


[Tx)  '         ^      ^    [TxJ 


The  asymptotic  behaviour  of  h!:   (x),  H^   (x)  changes  (Stokes'  phenomenon) 

A,         A, 

across  certain  branch  lines,  shown  as  thick  lines  in  Fig.  I5.   For  H   (x),  we 


ha 


ve  the  curves  h  (Re  a  =  0,  UmX  >   O)  and  h_  (Re(a  -  ±nX)    =  0,    JmX  <   O). 
These  ciorves  are  symmetrical  with  respect  to  the  origin  and  the  zeros  of  H   (x) 
are  asymptotically  located  on  them.   The  curve  h  cuts  the  real  axis  at  X,  =  x 
at  an  angle  of  jr/3.   The  tangent  to  this  curve  tends  to  the  vertical  direction 
for  I  x|  ->  00.   Asymptotically,  the  curve  approaches  X,  =  a|x|,  ti  ^  -it/2,  where  a 
and  T\   are  defined  by  (2.26)  and  (2.27)  with  p  replaced  by  x: 

a     -  exp[i(|  +  e)],      ti   =  €jU\2X/ex\.  (A. 5) 


o 

ll 
CJ 
<I> 

a: 


o 

V 

o 
a: 


o 

o 

1 

1 

o 

< 

< 

<   <M 

1 

II 

II 

II 

C'< 

ti^ 

^ 

X 

I 

-5 

o 

t 


O 

II 

o 
or 


H 
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(2 ) 
For  H^   (x),  we  have  the  branch  lines  h^(Re  a  =  0,  Cf-m  A,  <  0)  and 

A,  ^ 

h  (Re(a,  +  ±n\)   =  0,  CTm  A,  >  0 j  ,  which  are  complex  conjugate  to  h  and  h  , 

(2) 
respectively,  and  where  the  zeros  of  Yl        (x)  are  asymptotically  located. 

A. 

In  addition,  we  have  the  portions  of  the  real  axis  denoted  in  Fig.  15  by 

j   (from  -00  to  -x)  and  j  (from  -x  to  x),  where  the  zeros  of  J^(x)  are  located. 

These  curves  divide  the  \   plane  into  5  regions,  A  to  E  in  Fig.  15, 

and  the  asytnptotic  behaviour  of  h!;   (x),  h!;   (x),  J,  (x)  and  J  ^(x)  in  these 

A,         A.         A»  —  A# 

regions  is  shown  in  this  Figure.   Rote  that 

J  (x)   ->  (2rt^)"2(ex/2A)^       {\X\    ^oo)  (A.6) 

in  all  regions. 

For  each  function,  there  is  a  domain  where  it  tends  to  zero  for 
I  A.|  ^  «>,  whereas  it  tends  to  infinity  outside  of  this  domain.   For  J  (x),  this 

domain  is  region  A.   For  H,   (x),  it  is  the  domain  between  hp  and  the  cirrve 

(2 ) 
A,  =  -a|A,|,   1]  -^  -3Jt/2  in  the  lower  half -plane.   For  H^^   (x),  it  is  the  domain 

between  h-,  and  the  curve  A.  =  o|a.|,  t]  ->  3it/2  in  the  upper  half -plane.    Finally, 

J  (x)  -^  0  in  regions  C  and  D. 
—  A, 

These  resiilts  have  to  be  modified  in  the  neighbourhood  of  each  of  the 

branch  lines,  where  the  two  representations  for  the  same  function  on  different 

sides  become  of  comparable  order  of  magnitude.   We  then  must  take  for  the 

function  the  sum  of  the  two  representations.   This  is  indicated  by  the  shaded 

regions  in  Fig.  15. 

Thus,  we  have 
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H^^^x)  «  2Ae^"/^  sinh(a  -  i^)   in  AB ,  (A.  7) 

hJ^"'"^(x)  X  -2A  exp(-lnX,  -  if)sinh(a  -  inX  +  Ir)  in  DE,  (A. 8) 

H^^^(x)  «  -2A  exp(inA,  +  iv)sinh(a  -  ±nX  -   iv)   in  BC,  (A.  9) 

H^^^(x)  ;=;  2Ae"^"/  sinh(a  +  i|)   in  EA.  (A. 10) 


What  is  the  width  of  the  shaded  regions?   In  AB,  for  instance^  we  have 
H^   (x)  ^  A(e   -  ie   ),  while  outside  of  AB  one  of  the  two  terms  dominates, 
so  that  exp(2|Re  a| )  »  1.   Thus,  the  boundary  curves  of  AB,  shown  by  broken 
lines  in  Fig.  15,  can  be  defined  by 

Re  a   =  +  C,  (A.  11) 

where  C  is  a  constant  such  that 

e  "^  »  1,  (A.  12) 

-2C 
i.e.  e    may  be  neglected  within  the  required  degree  of  approximation. 

Asymptotically,  with  X.  =  o|\|,  we  find  that  (A. 11)  corresponds  to 

(^  +  ^1)1^1  =  +  c  ,  (A.13) 

where  a   and  r\   are  given  by  (A.  5).   Thus,  the  angular  width  of  the  region  AB  is 

Ag   =  2C/\x\  A,i2\x\/ex)  (A.li^) 

and  the  corresponding  arc  length  |X.|Ae  tends  to  zero  like  (/»i/|a.|  )   .   Similar 
results  are  valid  for  the  other  shaded  regions. 
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It  must  be  noted  that,  due  to  our  choice  of  the  phase  in  {A. 3),   A  goes 
over  into  -A  and  the  phase  of  a,  changes  by  2±nX   on  crossing  the  line  j',  so 

that,  in  spite  of  appearances  to  the  contrary,  the  representations  for  h!:   (x), 

(2) 
E\      (x)  and  J  ,(x)  given  in  Fig.  15  are  continuous  across  j',  while  that  for 

J  (x)  changes.   In  particular,  on  j'  itself,  with  \   =  -n,  we  find 
A, 

J_  (x)  X     A(|a,x)  <  sin(nM)exp[-a(|a,x)J  +  ^  cos(7tfi  )exp[a(n,x )] 

(n  >   x),       (A. 15) 

where  A  and  a,  are  given  by  (A.l)  and  (A. 2).   Thus,  J,(x)  has  infinitely  many 
zeros  on  j ' ,  located  asymptotically  very  close  to  the  negative  integers. 

The  asymptotic  expansions  given  in  Fig.  15  should  be  employed  for 
I  A,|  »  X.   For  smaller  values  of  \x\,    additional  regions  have  to  be  considered 
(25).   We  shall  require  only  a  few  additional  results. 

In  region  BE,  in  the  neighbourhood  of  the  real  axis,  we  may  employ  the 

Debye  asymptotic  expansions 

1 
(x^  -    X'~  )^   -    X   cos  "  - 


H^^'2^(x)   =   (2/«)2(x2  -  ;.^)"^/^xp<±  i 


2   ,2,t  _  ^  ^^^-1  X  _  n 

X    k 


1  + 


-,  fi^^^^.-.V- 


/.  ,2)4  V    3  (,2.  ^2)     y 


,   (A. 16) 


where  the  upper  signs  refer  to  hI;   (x)  and  the  lower  ones  to  e\      (x),  and 

A,  A, 

(x^  -  X,^)""^/   ■>  0,   0  <  cos"-'"(Vx)  <  Jt/2  for   -x  <  A.  <  x. 

These  expansions  fail  in  the  neighbourhood  ofX,=  +x.   If|A,-x| 

1/? 
becomes  comparable  with  |  A.|  '  ,  we  must  employ  the  expansions  (3_^,  pp.  367?  hkS) 
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h[^'^^(x)   =  2  exp(+ijt/3)(2A)^/3^i[exp(±2in/3)(2A)^/3^X-x)]  +  Cf'(\"^),   (A.l?) 
J  (x)   =  {2/x)^/^  A±l(2/x)^^^iX-^)']    +    ^ix'^),  (A.18) 


where  Ai(z)  denotes  the  Airy  function,  defined  in  Appendix  D. 
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Appendix  B.   Asymptotic  behavioiir  of  g(?c,P,p) 
In  order  to  prove  (U.2),  we  note  first  that,  according  to  (2.9)  and  (2.15), 

g(-^,P,p)   =  g(X.,P,P),  (B.l) 

so  that  it  suffices  to  study  the  asymptotic  behaviour  in  the  upper  half -plane. 
For  this  purpose,  we  shall  subdivide  it  into  two  overlapping  regions,  as  shown 
in  Fig.  16. 

In  region  1,  we  may  rewrite  g(>i., 3,p)  as 

g(^,P,p)   =  2[j^(P)  H^2)(p)  _  E^^\&)    J^(p)]  .  (B.2) 

(2) 
Substituting  J  and  H^   by  their  asymptotic  expansions,  valid  in  regions  A, 

AB  and  B  of  Fig.  15,  we  find 

g(X.,P,p)  :=  iA(X,P)A(X,p)  <^exp[a(\,p)  -  a(\,p)J  - 

-  exp[-a(X,p)  +a(X.,p)]>  .  (B.3) 

In  region  2,  g(A,, P,p)  may  be  rewritten  as 

g(\,P,p)   =  2e^"^[H[^^(P)  J_^(p)  -  J_^(P)  H^^^(p)]  .       (B.l^) 

Substituting  IT,   and  J_,  by  their  asymptotic  expansions,  valid  in  regions  C, 
BC  and  B  of  Fig.  15,  we  are  led  again  to  the  same  result  (B.3).   Taking  into 
account  (A. U),  this  yields 

(B.5) 


«(^.M)  ~~    i[(f)    -(f) 


ImX 


ReX 


Fig.    16 
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for  I  X.|  ->  00  in  the  whole  upper  half -plane.  It  then  follows  from  (B.l)  that 
the  same  resiilt  is  valid  in  the  lower  half -plane.  This  proves  the  validity 
of  (^.2)  in  the  whole  X   plane. 
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Appendix  C.   Asymptotic  behavioior  of  the  Legendre  fimctions 


The  functions  Q   (cos  0)  and  Q    (cos  0)  employed  in  the  present  work 


are  defined  by 


.(1,2) 


2i 


(cos  9)   =  I  [P  (cos  9)  +  —  Q  (cos  9)],    (C.l 


where  P  and  Q  are  Legendre  functions  of  the  first  and  second  kind,  respectively. 
We  have  the  following  relations  (35): 


P   (cos   9 


Q;   i(-cos   9 

A.- 2 


iHa 


Q;    l(-cos   0 
A,- 2 


P^   i(-cos    0 

A,- 2 


P^     i(-COS     0 
A,- 2 


Qi^'icos   9)   +  q'-^''(cos   9), 


=     le  '^\    L  cos    9), 

A,- 2 


.    ln\  ^( 1) /  „N 

-le  Q:    i(cos   9), 


(C.2) 

(C.3) 
(C.4) 


ig-i«^  P     j^(cos    9)    -  2i   cos(n\)Q;"^i(cos   9),      (C.5) 


-ie^"^  P,    i(cos   0)   +  2i   cos(rtX,)Qrl(cos   9).        (C.6) 

A,"  2  A,"  2 


Both  q!;  1  and  q!;  i   have  poles  at  the  negative  half -integers,  which  are  cancelled 

A,- 2         "-~2 

in  P  \. 

A."2 

If  e  <  9  <  jt  -e,  I  X,|  »  1,  |x,|€  »  1,  the  following  asymptotic  expansions 
are  valid  (35,  PP-  237,  2l+0): 


Q$V^(cos   9)      = 

A,- 2 


(2rtX  sin  9)2 


P^   i(cos   9) 

A.- 2 


1  ±  ^-§f-^  +  e-o.-^: 

rt\      ,    cot    0 


(c.t; 


(7rk^)'h(^^-!)  ^^  sin(.0-f) 


^       sin  f\0  -  |]  +    er(x"^) 


(c. 


97  - 


If  0  <  0  <  e,  |a,|£  <  1,  I  A-l  »  1,  we  may  employ  the  expansion  (35,  V-    ^^3) 


2  0 
P^  i_(cos  0)   =  J  (u)  +  sin  - 

A."  2  o  ^ 


u 


J^Cu)  1 


6  -^3^^^  "  ^2^^^  ^   2u 


+  &    (sin^l  )  ,    (C.9: 


where 


u  =  2A,  sm  — 


(CIO) 


.   8 


A  uniform  asymptotic  expansion  of  the  Legendre  function  has  been  given 
by  Szego  (36): 


P,  i(cos  9)   =  (-r^Q  ' 
X-2  \sin  0 


J  (X9)         - 

j^(x.0)  +  ^(0  cot  0- 1)  -^^ —  +  e'  (x    ) 


;c.ii) 


This  reduces  to  (C.8)  when  |x|0  »  1  and  remains  valid  for  0  -*  0. 
In  Section  IX  we  made  use  of  the  inequality 

|P.  _i(cos  0)|   <  (cos  9)  2exp(rt|i/2)        (^  >  0,      0  <  0  <  7t/2).  (C.12; 

To  prove  this  inequality,  we  start  from  the  integral  representation  (25,  p.  38?) 

CO 

1 


P  (cos  0) 


r(v+i: 


exp(-t  cos  0)  J  (t  sin  0)t  dt 


(0  <  0  <  n/2,   Re(v+1)  >  O), 

(C.13) 


where   r(z)   is  the  Gamma   function.      It   follows  that 

00 

1        r 


Ip.      i(cos   0) 

11-1-2 


|r(i^+i)|    J 


exp(-t   cos   9)   t      dt 


(cos   9)   2 


xCil 


'r(in+i)| 


cos^h   (jtm)"I2  e"^^^ 


cos    0 


(cos   9] 


T     (0  <  9  <^)   , 


which  proves  (C.12, 


Appendix  D.   The  Airy  Function 


The  Airy  function  is  defined  by 

.1/3 


Ai(z)   =  ^^  j     cos(t3  +  3^'^^zt)dt  .  (D.l) 

o 

It  can  also  be  expressed  in  terms  of  Bessel  functions  of  order  I/3  (3^1,  p.  hk6] 


We  have 


¥[Ai(z),   Ai(ze-2i«/3)]   =  e^^"/*^/2«,  (D.2) 

where  W  denotes  the  Wronskian.   Also, 

Ai(z)  +  e2i«/3  AKze^i'^/S)  -.  e'^i'^/S  j,,^,,-2±n/3^      =  0.     (D.3) 
The  asymptotic  expansion  of  Ai(z)  for  large  |z|  is  given  by  {3}^,   p.  ^^) 

Ai(z)  ;^    1  ./,  [1  +  0  (T^)]         (|argz|<Tt),      (D.^) 

Ai(-z)  ^  rt-2  z-l/^Jsin(^  +  f)[l  +  ©'(r^)]  -cos(C+f)  e^(r^)> 

(jarg  z|  <  2it/3),    (D.5) 

where 

^  =  I  z3/2  .  (D.6) 

The  zeros  of  the  Airy  function  are  all  located  on  the  negative  real  axis. 
If  -x  denotes  the  nth  zero,  we  have,  for  large  n  {3k,   p.  U50), 


^n  =   [f  (n  -  t)]'^^[l  +  »'(^'^)]        (n»l),   (D.7) 
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Ai'(-x  ) 

n 


(-ir'  .-* 


^  (n 


1, 
4' 


11/6 


(n  »  1)     (D.8) 


lYie   first  five  values  of  x  and  the  corresponding  values  of  Ai'(-x  )  are 
listed  in  Table  I  (3^,  p.  ^78). 
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Appendix  E.   The  Lommel  Functions 


Lommel's  functions  of  two  variables  (25,  pp.  537-550)  are  defined  by 

2m-v+2 


Vju,v)   =  cos(^-  .-  -.-J  + 


(-iT  (^) 


■^2m-v+2^^)-   ^E-i: 


m=0 


In  particular. 


V  (u,u)   =  p[j  (u)  +  cos  u]. 


V-j_(-u,u)   =  -  -  sin  u, 


Vq(u,0)   =  1, 


V^(u,0)   =  0 


:e.2) 


(E.3) 


The  following  integral  representation  is  valid; 

00 

v-1    P 


V  (u,v)  +  iV  ^ (u,v 
V  v-1  ' 


v+2. 


v-2 


J   (vt)exp[iu(l-  t  )/2]  t  ^  dt 


u 


(u  >  0,  V  >  0,  Re  V  >  -).  (E.4) 


In  particular. 


J^(vt)exp(iut^/2)tdt   =   (i/u)exp(iu/2)[v^(u,v)  +  iV^(u,v)]  .     (E.5) 


For  large  |u[  and  fixed  v  and  v,  we  may  employ  the  asymptotic  expansion 

{E.6) 


V^(.,v)  .   ^(-ir(v/u)->^.^(v). 


m=0 
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Appendix  F.   Evaluation  of  f  (k^jt  -e) 


r 


To  eval-uate  the  first  integral  in  (9.51)^  we  may  employ  the  expansion 

oo 

tan(jrax)      =     i   +  2i   y      (-1)'^  exp(2injta  x),  (F.l) 


n=l 

which  gives 


00  <» 

/     exp(-x    )J   (tox)tan(nax)xdx     =•  i    /   exp(-x    )J  (u)x)xdx 


+  2i   y      (-1)"      /  exp(-x     +  2inrtax)J^(oJx)xdx.  (F.2) 

n=l  o 

According  to  Weber's   integral  formula   (25,   p.    393), 

CO 

/2  12 

exp(-x    )J   (u)x)xdx     =     -  exp(-oj  /k).  (F.3) 

o 

On  the   other  hand,   hy  partial   integration, 

00 

r  ?  1-2 

/     exp(-x     +  2innax)J   (cox)xdx     =     -  5     +     ©^(P      ),        (F.U) 

J  °  (2njra)^ 

o 

so  that,  finally, 

00 

/  exp(-x^)j  (LJx)tan(rtax)xdx  =  |  exp(-w^A)  +  -^  +  O'Cp"^),   (F.5) 
o         °  21^ 

where  we  have  employed  the  well-known  formula 

,   >n+l        2 

^-^^   =   15.  (F.6) 

n=l 


102 


It  is  clear  from  the  above  calculation  that,  since  we  are  neglecting 

_2 
terras  of  the  order  of  p   ,  we  may  replace  tan(rtax)  by  i  in  all  remaining 

integrals  of  (9.51)-   To  compute  these  integrals,  we  employ  Hankel's  formula 

(25,  p.  393): 


,   2 .  ^  ,   .  m-1 ,       n  o  ^ 
exp(-x  )J^(cjxjx   dx  =   — 


f^  iff   exp(-u.2/4)  0  (^  .  1,  n.  1,  ^j)  ,   (F.7 


where  0(a,b,z)  is  Kummer's  confluent  hypergeometric  function. 


Substituting  the  above  results  in  (9.51),  we  get 


f^(k,"-0 


exp(-w  /h)   +  


2 
12  a -J 


+  ^   exp(-c.7^)$(-2,l,^; 


,1.26     ,   2/,  . 
+  |j;sin  -  exp(-w  /h) 


r  h  2     2 

^  0(1,^,—  )  -  — *(l,3,— )  +  ^(1,2,—  . 


2 

1+ 


2 
If 


+  &'(P"^) 


(F.8; 


It  follows  from  the  definition  of  ^(a^bjz)  that 


^(-2,1,—  ) 


2    k 

OJ      CO 

2  ^  32  • 


(F.9) 


On  the  other  hand,  we  have  (38) 


<l'(l,n  +  l,z)   =  ne  z  r(n,z), 


(F.IO) 


where  Y"(n,z)  is  the  incomplete  gamma  function,  and 


n-1 


r(n,z)   =  (n-1):   1  -  e 


-z 


m=0 


(F.ll) 
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so  that 


n-1 


$(l,n  +  l,z)   =  ^  f  e^ 


m 


(F.12) 


m=0 


Substituting  (F.9)  and  (F.12)  in  (F.8)  and  taking  into  account  (9.52),  we 
finally  get 


f^(k,it-  < 


J  exp  (-2ip 


.  k 


i-i-g?*«'(^   ) 


:f.i3) 


-  lOif  - 


Footnotes 

This  transformation  seems  to  have  been  first  employed  by  Bremmer 
(2,   p.  210).   It  was  subsequently  adopted  by  several  authors. 

2 

An  integral  representation  related  to  (2.3^)^  but  involving  trigonometric 

instead  of  Legendre  functions,  was  given  by  Franz  and  Galle  (8). 

For  a  discussion  of  the  width  of  this  neighbourhood,  see  Appendix  A. 

k 

Actually,  the  expression  within  square  brackets  in  the  second  term 

of  (6.23)  differs  from  Keller,  Lewis  and  Seckler's  by  the  powers  of  2  in  the 

denominators,  which  are  all  less  by  one  unit.   However,  it  can  be  verified 

that  this  is  due  to  a  misprint  in  their  paper. 

^  Cf.  Eqs.  (26)  -  (27)  and  the  related  discussion  on  p.  328  of  (3I)  and 
Eqs.  (5)  -  (8)  and  {k8)   of  (32). 

Similar  figures  appear  in  (33  )•   However,  they  contain  several  mistakes. 

Franz  (5,  p.  36)  incorrectly  states  that  E  and  B  are  the  domains  where 
H;"'"^(x)  ^0,  e[^\x)   ^0,  respectively. 

A*  A. 

Q 

Cf.  also  (37).   However,  there  are  several  mistakes  in  the  expansion 
given  in  this  reference. 
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Table  I.   The  First  Five  Zeros  of  Ai(-x) 


n 

X 

n 

Ai'(-x 
^  n 

1 

2.33811 

+0.70121 

2 

U. 08795 

-0.80311 

3 

5.52056 

+0.86520 

h 

6.78671 

-0. 91085 

5 

7. 9^i^l3 

+0.9473^ 

Flgtire  Legends 

Fig.  1.   Paths  of  integration  in  the  X  plane. 

Fig.  2.   The  path  of  integration  in  (2.19)  must  begin  and  end  at  infinity  to 
the  right  of  the  shaded  regions.   A  possible  path  that  is  symmetrical  about 
the  origin  is  shown  (0  <  n/2). 

Fig.  3.   The  path  of  integration  in  (2.20)  must  begin  and  end  at  infinity  to 
the  right  of  the  shaded  regions.   A  pos'sible  path  that  is  symmetrical  about 
the  origin  is  shown  (0  >  Jt/2). 

Fig.  h.      XXX  Poles  of  the  S-function,  located  on  curve  h  .   The  contour 
C  passes  half-way  between  consecutive  poles. 

Fig.  5.   Diffracted  rays  T  T ' P  and  TpT'P  reaching  a  point  P  in  the  geometrical 
shadow  region. 

Fig.  6.  XXX  Poles  of  the  S-function;  o  o  o  zeros  of  the  S-function;  •  •  • 
poles  of  Q^  X(cos  0).   The  path  of  integration  T   goes  through  the  saddle 

A,- 2 

points  A.  =  kp  and  X  =  psin0   (0<  3t/2). 

Fig.  7.   Physical  interpretation  of  the  left  saddle  point:  p  is  the  impact 
parameter  of  the  incident  ray  AB  which  reaches  point  P  after  geometrical 
reflection  at  the  surface. 

Fig.  8.   Diffracted  rays  T  T' P  and  TgT'P  reaching  a  point  P  in  the  lit  region. 

Fig.  9.   I^ths  of  integration  for  (7.13);  w  is  the  saddle  point,   x  x  x  poles 
of  rAi(w)]"  .   The  path  must  begin  and  end  at  infinity  outside  of  the  shaded 
sector. 


Figure  Legends,  continued 

Fig.  10.   Path  of  integration  in  (8.23).   x  x  x  Poles  of  the  integrand.   The 
integrand  goes  to  zero  at  infinity  outside  of  the  shaded  sector,  except  along 
the  line  of  poles. 

Fig.  11.   Geometrically  reflected  ray  in  direction  9. 

Fig.  12.   Ifeths  of  integration  for  (9.8)-   x  x  x  Poles  of  S(A.,^). 

Fig.  13.   Paths  of  integration  in  the  X  plane,   x  x  x  Poles  of  S(A., P) 

Fig.  1^.   Division  into  regions  (the  angles  are  greatly  exaggerated).   The 
separation  between  regions  is  indicated  by  the  broken  lines. 

Fig.  15.   Asymptotic  behaviour  of  the  cylindrical  functions  Z  (x)  in  the  index 

K 

plane  (x  >  0).   The  zeros  of  H^   (x)  are  asymptotically  located  on  curves  h 

A,  J_ 

(2) 

and  h  ;  those  of  H^   (x),  on  hp  and  h_p;  those  of  J  (x),  on  j  and  j'. 


Fig.  16.   Division  of  the  X,  plane  into  regions  ^s^\^  Region  1   V////    Region  2 


The  two  regions  overlap  in  the  cross-hatched  domain.   The  curves  in  broken 
line  correspond  to  those  of  Fig.  15- 


List  of  symbols  (in  order  of  appearance) 


Tp  =  I.e.  greek  "psi" 

e  =  1. c.  greek  "theta" 

p  =  I.e.  greek  "rho" 

P  =  I.e.  greek  "beta" 

X  =  1.  c.  greek  "lambda" 

jt  =1.0.  greek  "pi" 

e  =  I.e.  greek  "epsilon" 

cr  =  I.e.  greek  "sigma" 

T)  =  1.  c.  greet  "eta" 

6  =  eap.  letter  "0"  (order  symbol) 

7  =  1. c.  greek  "gamma" 
5  =  I.e.  greek  "delta" 
r  -  cap.  greek  "gamma" 
a  =  1. e.  greek  "alpha" 
^  =  1. c.  greek  "zeta" 
X  =  greek  "chi" 

<jr  =  cap.  script  "ef" 

T  =  I.e.  greek  "tau" 

V  =  I.e.  greek  "nu" 

^  =  1. c.  greek  "csi" 

cp  =  1.  c.  greek  "phi" 

f  =  "ef  tilde" 

A  =  cap.  greek  "delta" 

(i  =  I.e.  greek  "mu" 

(^  =  1.  c.  greek  "omega" 

0  =  cap.  greek  "phi" 
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